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Samenvatting, 


In dit verslag bestuderen we de interactie van een accre= 
tieschijf met het magnetische dipoolmoment van een neutronen= 
ster, We gaan daarvoor eerst uitvoerig in op de methoden om de 
vergelijkingen voor de accretieschijf op te lossen en behande 
len vervolgens het model van Ghosh en Lamb, Ondanks dat dit moe 
del vrij algemeen geaccepteerd lijkt, moeten we toch een be- 
langrijke inconsistentie constateren, Ook de afleiding van het 
effectieve geleidingsvermogen is in dít model nogal dubieus, 

We hebben dan ook geprobeerd hiervoor een betere af leiding 
te vinden, Zo vonden we het accretieschijfmodel van Ichimaru, 
waar het effectieve geleidingsvermogen eenvoudig volgt uit de 
magnetische turbulentie die ook de anomale viscositeit bepaalt, 
Nieuwe berekeningen met dit effectieve geleidingsvermogen in 
het kader van het model van Ghosh en Lamb geven weliswaar quan= 
titatief andere resultaten, maar bevestigen opnieuw het bestaan 
van een evenwichtsperiode voor de neutronenster, waarop het to= 
tale koppel op de neutranenster nul wordt, Helaas blijken ook 
deze resultaten inconsistent en is er nog verder werk nodig om 
tot betrouwbare quantitatieve uitspraken opgrond van het model 
van Ghosh en Lamb te komen, We constateren, dat een eventuele 
relatie tussen het effectieve geleidingsvermogen en de anomale 
viscositeit hierbij speciale aandacht verdient, 

Vervolgens gaan we in op hét model van Anzer en Börner, 
waarin het magnetische dipoolmoment niet parallel ís aan, maar 
loodrecht op de de rotatieas van de neutronenster (en de accre= 
tieschijf) staat, Ook ín dít model constateren we een aantal 
onvolkomenheden, die nadere bestudering verdienen, 

Als laatste gaan we dan ín op het model van de scheve To= 
tator, waarin een oneindig goed geleidehde dunne schijf wordt 
bestudeerd in een zo algemeen mogelijke geometrie: de hoeken 
tussen alle twee de rotatieassen en het magnetische dipoolmo= 
ment zijn volledig willekeurig, Dit model blijkt naast het 
grillige periode gedrag onder andere ook de precessie van de 
accretieschijf (Her X=1) en van de neutronenster (55 433) te 
kunnen verklaren, Bovendien veranderen de resultaten nauwelijks, 
als de schijf niet ideaal geleidend ís maar een klein deel van 
de magnetische flux in de schijf kan doordringen, zoals in het 
model van Ghosh en Lamb, Dit model van de scheve rotator ver= 


dient voorlopig de voorkeur, hoewel nog veel onzeker blijft, 
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Ï Inleiding. 


Sinds de opkomst van de Röntgensterrenkunde in het begin 
van de zeventiger Jaren is een groot aantal Röntgenbronnen ín 
ons eigen melkwegstelsel ontdekt, Omdat een aantal van deze 
Röntgenbronnen geïdentificeerd kon worden met een eclipserende 
variabele of met een spectroscopische dubbelster, ontstond het 
idee, dat deze Röntgenbronnen het gevolg zijn van de massaover= 
dracht ín deze dubbelsterren van een begeleider naar een come 
pact object, De hoge Röntgenlichtkracht van deze bronnen, van 
de orde van 10° erg/sec, en de hoge temperaturen, van de orde 
van 108 K als de straling van thermische oorsprong is, wijzen 
ook op massaoverdracht naar een compact object, Het compacte 
object kan een witte dwerg, een neutronenster of een zwart gat 
zijn. De massaoverdracht kan het gevolg zijn van een sterrewind 
van de begeleider of van Roche=lob=overstroming ven de begeleie= 
der, als deze zijn binnenste kritische oppervlak vult (zie fie 
uur 1), Een goed overzicht van deze Röntgenbronnen wordt gege= 


ven in Lewin & van den Heuvel, 1985, 


a k | 
& Î JZ wind 








Figuur 1 Massaoverdracht in dubbelsterren, 
a) Accretie door de wind van de begeleider, 
b) Accretie door Roche=lob=overstroming. 
In het geval van Roche=lob=overstroming is het impulsmo= 
ment van de overstromende materie te groot om direct op het 


compacte object te vallen en zal zich een accretieschijf vormen, 





Van ongeveer 20 Röntgenbronnen is bekend, dat zij regel— 
matige pulsaties vertonen met periodes van één tot enkele hon= 
derden seconden dit zijn de Röntgenpulsars. Aangenomen wordt, 
dat het compacte object in een Röntgenpulsar een magnetische 
neutronenster (of een witte dwerg) is, waarvan het magnetische 
dipoolmoment niet parallel is aan het impulsmoment van de neum= 
tronenster (of de witte dwerg); als dit dipoolmoment groot ge= 
noeg is zal de accretie langs de magnetische veldlijnen en dus 
op de polen plaatsvinden, Voor een waarnemer in de juiste ge= 
zichtslijn zal een pool één maal per rotatie periode zichtbaar 
zijn; de waarnemer ziet dan dus Röntgenpulsen met een periade, 
die gelijk ís aan de rotatie periode van de neutronenster (zie 
figuur 2), Voorwaarde hiervoor is, dat nog voor de materie op 
de neutronenster (of de witte dwerg) valt, de materie-druk ge= 
lijk wordt aan de magnetische druk, De radiële afstand, waarop 
dit gebeurt, wordt de Alfvón=straal, Eys of ook wel de magne= 


tosferische straal, En? genoemd, Binnen deze straal is de ac= 


cretiestroming gekoppeld aan het magnetische veld, 


waarnemer 
vg 





Figuur 2 De vorming van Röntgenpulsen door accretie op de 
magnetische polen van een neutronenster (of een 
witte dwerg). 

De magnetische druk is sterk afhankelijk van de radiële 
afstand (m rx"), waardoor de preciese waarde van de Alfvén- 
straal weinig afhankelijk ís van de vorm van de accretie=stroom 


en typisch is 5 ‚ 
Ta= „oft Cacm)=1/7 n2/7 ë (1) 


Hierin is gy het magnetische dipoolmoment, G de gravitatie con= 
stante, M de massa van de neutronenster (of de witte dwerg) en 


M is de massaoverdracht per seconde, (zie Henrichs, 1983) 





Van de meeste Röntgenpulsars blijkt de periode te veran= 
deren. In figuur 3 ís voor acht Röntgenpulsars de periode uit= 
gezet tegen de tijd, Duidelijk is dat vier van de acht een re= 
gelmatige afname van de periode ('spinup!) vertonen, Het peri= 


ode=gedrag van de andere vier ís echter erg onregelmatig, 


Pulse Period (sec} 








3 : 
Fig. (Adapted from Rappaport & Joss 1982) Pulse period histories for 
8 X-ray pulsars. In cases of a more or less regular trend a dashed line 


is fitted to the data points. 


Een regelmatige af name van de periode kan als volgt ver= 
klaard worden (zie Rappaport & Joss, 1977): 

De materie die naar de neutronenster (of de witte dwerg) 
overstroomt, bezit impulsmoment en vanaf het moment, dat de ma= 
teriestroming aan het magnetische veld ís gekoppeld, ís het to= 
tale impulsmoment van de materie en de neutronenster behouden, 
In het geval van een accretieschijf betekent dit een impulsmo= 
mentoverdracht van 


3 = (emr‚)£ m Ee, (2) 


Voor de neutronenster geldt dean een frequentie toename van 











° | e 

Q= (enr)E m1, (3) 
Hierbij ís verondersteld, dat het traagheidsmoment, I, van de 
neutronenster constant is, Deze aanname is goed, als de massa 
toename van de neutronenster kléín is ten op zichte van zijn 
massas voor vrijwel alle Röntgenpulsars is dit in orde, Als we 
veronderstellen, dat alle gravitatieënergie van de invallende 


materie wordt omgezet in Röntgenstraling, vinden we de volgende 


relatie tussen M en de Röntgenlichtkracht, Le 
Gm À 
Loe Me (4) 
Uit vergelijkingen (3) en (4) vinden we een relatie tussen de 


periode af name, =P, en de Röntgenlichtkracht: 


els vet 2/7 sf efr 6/7 
= (0,7 107* 3171) ie „217 „sd R/ L$/ p (5) 


1 
Du. 


Hierin is I het traagheidsmoment in Tin or.en°, Man het mag= 


A9 30 3 
netische dipoolmoment in 10 Gem 


‚ m de massa van de neutro= 
nenster in 1,áé Mg Re de straal van de neutronenster in 106 cm 
en is ks, de Röntgenlichtkracht in 1057 ero/sec. 

In figuur 4 staan de -P/P's uitgezet tegen de pLS/?1s van 


een aantal Röntgenpulsars, De door=getrokken lijn in deze fie 
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Fig.3.2. (Adapted from Rappaport & Joss 1977) According to & simple 
accretion-torque model the time-scale for the change in rotation period 
(units yr) is proportional to P Le77 (see eq. 3.9.7). A best fit: 
through the data points (where Her X-1 and the two limits were excluded) ' 
is indicated by the solid line, The derived proporttonality factor 
agrees with the value as expected for neutron stars, whereas this factor | 
is expected to be much lower ' for. white dwarfs (dashed line). This, 
agreement gives very strong evidence in favour of the interpretation 


that X-ray pulsars are neutron stars, rather than white dwarfs. 
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guur is de beste aanpassing door de punten ( Her Xe1 en de twee 
limieten werden niet meegenomen), De hieruit afgeleidde evenre= 
digheidsconstante komt goed overeen, met wat we voor een neutro= 
nenster zouden verwachten, De onderbroken lijn geeft de ver= 

‘ wachte relatie weersals het compacte object een witte dwerg zou 
zijn, Op grond hiervan kunnen we concluderen, dat Röntgenpul= 
sars inderdaad magnetische neutronenesterren zijn en geen witte 
dwergen. Opvallend is echter de uitzonderlijke positie van Her 
X=1; de periode=af name is veel geringer dan de door de theorie 
‘voorspelde waarde, 

In de bovenstaande afleiding ís in feite verondersteld, 
dat de grenslaag tussen de accretieschijf en de magnetosfeer 
oneindig dun is. In werkelijkheid moet deze laag een eindige 
dikte, 6, hebben, waarbinnen de rotatie wordt afgeremd van Kepe 
Ierse rotatie aan de buitenrand tot corotatie met de magneto= 
sfeer aan de binnenrand van de grenslaag, In Ghosh & Lamb, 1979 
a en b, (in het vervolg G&L 1 resp, II) werken de auteurs een 
model uit, waarin hiermee rekening wordt gehouden, Twee belang= 
rijke grootheden in hun model zijn de corotatiestraal, dat is 
de radiële afstand tot de neutronenster, waar de Keplerse rota= 
tieperiode in de accretieschijf gelijk is aan de rotatieperiode 
van de neutronenster (N.8, de magnetosfeer wordt verondersteld 
een starre rotator te zijn), 

aa Sk sl 1/3 En | EL 
NS 
en de 'fastnessparameter', dat is het quotiënt van de rotatie= 
periode van de neutronenster en de Keplerse rotatieperiode op 
de magnetosferische straal (nu de straal van de binnenrand van 
de grenslaag tussen de schijf en de magnetosfeer), 
Ss 
a lr) « (7) 


« 


De dikte van de grenslaag blijkt in het model van Ghosh en 
„Lamb klein te zijn, deweZe Sar. Bovendien is de grenslaag 
‚niet in staat om de rest van de accretieschijf volledig van de 
magnetische veldlijnen af te schermen, De magnetische veldlijn= 
en penetreren over een uitgestrekt overgangsgebied, dat zich 
‘tot de radiële afstand Te uitstrekt, in de accretieschijf (zie 


figuur 5). 








„—MAGNETOSPHERIC TRANSITION UNPERTURBED_____ 
glen DISK F! LOW 
TRANSITION 


OO ATER 


Fig. Bi. (adapted from Ghosh & Lamb 1979b) Schematic picture of a 
possible configuration of an accretion disk around a magnetic neutron 
„ Star, In this picture the magnetic field Ís not completely screened by 
currents in the magretospheric boundary and hence the field lines may 
thread the disk. This leads to an additional contribution to the 
accretíon torque (see text and fig. Je i 


In dit model is het koppel op de neutronenster als gevolg 
van de grenslaag gelijk aan het koppel, dat uit de hiervoor ge= 
geven simpele afleiding volgt, omdat de visceuze spanning daar 
te verwaarlozen is en de grenslaag relatief dun is, Als de neu= 
tronenster langzaam roteert, deweze als WT 0, zal in het hele 
overgangsgebied de materie in de accretieschijf sneller roteren 
dan de magnetosfeers de magnetische veldlijnen worden dan mee 
gesleept, wat resulteert ín eenextra koppel op de neutronen 
ster, Bij een snel roterende neutronenster zal de corotatie- 
straal echter in het overgangsgebied kunnen liggen, De rotatie= 
snelheid als functie van de radiële afstand in zo'n geval is in 


fiouur 6 geschetst, 





Figuur 6 De rotatiesnelheid in de accretieschijf, 
Omdat de materie in gebid A sneller roteert dan de 
neutronenster, geeft dít een positieve bijdrage aan 
het totale koppels evenzo geeft B een negatieve bij- 
drage aan het totale koppel, Zie verder de tekst. 
(naar Henrichs, 1983) 
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In gebied A in figuur 6 toteert de materie in de accretieschijf 
sneller dan de magnetosfeer en resulteert de magnetische mee= 
sleping in een positieve bijdrage aan het totale koppel, In ge= 
bied B, tussen de corotatiestraal en de grens van het over= 
gangsgebied, Ps roteert de materie in de accretieschijf lang= 
zamer dan de magnetosfeer; hier resulteert de magnetische mee= 
sleping in een negatieve bijdrage aan het totale koppel, 

De vraag. is nu natuurlijk hoe groot de verschillende bij= 
dragen zijn en wat het totale koppel is, Om deze vraag te kun= 
nen beantwoorden is een quantitatieve oplossing van de stroming 
in de accretieschijf in het overgangsgebied en de grenslaag 
nodig, Deze wordt gegeven in G&L Te (zie ook Henrichs, 1983) 

De opzet van dit verslag is nu als volgt: 

We zullen eerst bespreken hoe de normale (zonder B-veld) bere= 
keningen voor de stroming in een accretieschijf worden gedaan, 
de vergelijkingen opstellen voor het geval er wel een magne= 
tisch veld is en de methode van oplossen van Ghosh en Lamb be= 
schrijven en hun resultaten geven, Daarna zullen we hun methode 
en‘hun resultaten kritisch bekijken en concluderen, dat een 
aantal aannames op zijn minst dubieus zijn en dat de resultaten 
inconsistent zijn. Vervolgens zullen we enige alternatieven 
aandragen en voor één nieuwe set aannames hun analyse gedeelte= 
lijk herhalen, Dan zullen we een discussie geven van de resul 
taten en twee alternatieve modellen uit de litteratuur kort be= 
spreken, Tot slot willen we proberen eventuele lijnen voor 


toekomstig onderzoek aan te gevene- 





II=1 Standaard accretieschijfmodellene 


Omdat een standaard schijfmodel de basis vormt van het 
model van Ghosh en Lamb, geven we eerst een kort overzicht van 
de theorie van accret{ieschijven, Omdat de problemen, die in de 
theorie over accretieschíijven nog bestaan, door werken in elke 
andere theorie, die de schijfmodellen als basis hebben, willen 
we op de bespreking daarvan de nadruk legger, Voor meer details 
willen we verwijzen naar Verbunt, 1982; we zullen onze bespre= 
king op deze uitstekende review baseren, 

Het is niet vanzelfsprekend, dat massaoverdracht ín een 
dubbelster leidt tot de vorming van een accretieschijf; hier= 
voor zijn twee noodzakelijke voorwaarden: 

1) de specifieke, kinetische energie moet klein genoeg zijn om 
gravitationele vangst door het compacte object mogelijk te 
maken, 

2) het specifieke impulsmoment moet groot genoeg zijn om direce 
te inval op het compacte object te voorkomen, 

Bij een magnetische neutronenster betekent de tweede voorwaarde 

dat het specifieke impulsmoment groter moet zijn dan het specie 

fieke impulsmoment van Keplers roterende materie op de magneto= 
sferische straal. 

Voor massaoverdracht door Roche=lob=overstroming zal E 
meestal aan de beide voorwaarden voldaan zijn en kan er zich 
dus een accretieschijf vormen. In figuur 7 îs de vorming van 
zo'n accretieschijf schematisch weergegeven, Omdat de drukgra= 
diënten ín de stroming vanuit het binnenste Lagrange=punt klein 
zijn, is de stroming vrijwel ballistisch (zie Lubow & Shu, 
1975); de materie valt niet direct op het compacte object, maar 
volgt een sterk excentrische baan rond het compacte object om 
weer tegen de inkomende stroom te botsen, figuur 7a, Daar vormt 
zich een schok, die de supersonische stroom stroomopwaarts niet 
kan beïnvloeden, maar stroomafwaarts de materie ín een cirkel- 
vormige ring doet;stromen, figuur 7b, Onder invloed van de vis= 
cositeit verbreedt de ring zich, waarbij de meeste materie naar 
binnen en een klein deel van de materie en impulsmoment naar 
buiten worden getransporteerd, Dit laatste vanwege het behoud 
van impulsmoment. (zie figuur 7c) Be expansie van de ring gaat 
door, totdat deze een straal dichtbij de Roehe-lob heeft bee 
reikt en de getijdewisselwerking met de begeleider het impuls= 


moment, dat naar buiten îs getransporteerd, over kan brengen 





a) initial gas 
stream 


b) formation of 
ring 


c) ring spreads 





d) disk is formed 


dî) side view 





Figuur 7 Een schematische weergave van de vorming van een 
accretieschijf (naar Verbunt, 1982). 
naar het impulsmoment van de baanbeweging (zie Papaloizou & 
Pringle, 1977).en totdat de binnenkant van de schijf de magne= 
tosferische straal heeft bereikt; de accretieschijf is gevormd, 
figuur 7d, Er is echter nog geen volledige duidelijkheid over 
de buitenste grens van de accretieschijf (zie beve Shu & Lubow, 
1981). We zullen hier echter verder niet op ingaän, omdat voor 
ons doel slechts het’ binnenste van de accretieschijf van belang 
is. 

Vanwege de hoge dichtheden en de relatief kleine vrije 
weglengte in een accretieschijf is een hydrodynamische be= 
schrijving van de accretieschijf noodzakelijk. Omdat het boven= 
dien de viscositeit ís, die de materie naar:binnen en het im- 
pulsmoment naar buiten transporteerd, moeten we de hydrodyna= 


mische vergelijkingen voor een visceuze vloeistof gebruiken, 
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Dat zijn: 
1) Behoud van massa dt Vv. (ev) (1) 
2) Bewegingsvergelijking 
e+ &.05) 


3) Energievergelijking 


zov.) - EV) - và (5) 


eE + (te) 

Hier in is $ de zwaartekrachtspotentiaal, À de visceuze span= 

ningstensor, € de specifieke inwendige energie en is À de - 

energief lux als gevolg van warmte geleiding en straling, 
Met de volgende aannames kunnen we deze vergelijkingen 

_ vereenvoudigen: 

1) Axisymmetrie « Sp = cai 
Dit betekent, dat de invloed van de begeleider wordt ver= 
waarloosd, Merk op, dat we cylindrische coördinaten r‚ & en 
z gebruiken, 

2) De accretieschijf is dun: h/r«&1. 

Hier is h de halve dikte van de aeccretieschijf, Dit houdt 
in, dat de druk-krachten in de schijf kleiner zijn dan de 
gravitatiekrachten, dan is p «eve. Bovendien volgt nu uit 
vergelijking (1), dat ve (h/r)v, en dus dat v‚«v_. 

3) De elementen van de visceuze spanningstensor zijn kleiner 
dan de druk: Er sP en dus t‚j cer. Dit betekent, dat de 
viscositeit leidt tot een kleine radiële snelheid, zie bev, 
vergelijking (7): vk Vge | 

4) De accretieschijf is stationair: de = Oy 

In het vervolg duiden we de dichtheid, druk en temperatuur 
in het centrale vlak, z = 0, aan met respectievelijk @‚ P en Te 


We definiëren de halve dikte van de accretieschijf, h, nu door: 
«a 


ho = (z) dz . (4) 
e=! ez) 
Integratie van vergelijking 1 over z en Ò geeft ru: 

à.: á 

Sr rehv‚, ) = 0, (5) 
De r}) component van vergelijking (2) wordt: 

v‚ = EE ‚ dei, Keplerse rotatie, (6) 


Evenzo de &) component van vergelijking (2): 


rehv, rl Dg ) = erf E tok 1 . (7) 





1 


waarin T ai s { tz) dz. In cylindrische coördinaten geldt 


voor de re-component van de visceuze spanningstensor (zie 


Landau & Lifshitz, En nd ane 15,15): 


eten teint. (8) 


Hierin is \ de dynamische viscositeitsconstante, 
Vervolgens geeft de z) component van vergelijking (2) de verge= 
lijking voor hydrostatisch evenwicht: 
2 Er, 

p = $@h (9) 
Tot slot vinden we uit de over z geïntegreerde energievergelij= 
king, dat de energieproductie als gevolg van de viscositeit, 
ri gelijk is aan de energief lux door het oppervlak van de 
schijf, Q': ee 

wt Öl ie tenn dt 

Q =| 5 dz Q= Tr re). (10) 
Op eenvoudige wijze kunnen uit de bovenstaande vergelijkingen 
de set 'standaardeschijf ‘evergelijkingen worden afgeleid (zie 


Verbunt, 1982): 


n 
rehv, = = zn (11) 
2 Gm 
Ve T (6) 
P, 2 
Mrvg (1 VE ) = 47r haa (12) 
p= hen” EE (9) 
8 
- + GMM Fy, 
ee (af). (13) 
anr 


De factor ( VvE ) is het gevolg van de keuze voor de rande= 
voorwaarden aan de binnengrens, T,, van de accretieschijf, Om= 
dat deze factor aanleiding geeft tot singulariteiten in de op= 
lossing en de juiste keuze voor de randvoorwaarden in feite 
onbekend is; zullen we deze factor in het vervolg weglaten, 

Om de set stándaardeschijf-vergelijkingen op te lossen 
zijn nog drie extra vergelijkingen nodig, Dat zijn een toe= 
standsvergelijking, een vergelijking voor Q” en een vergelij- 
king, die de dynamische viscositeitsconstante definiëert, Voor 
de toestandsvergelijking is de volgende keuze voor de hand 
liggend: ij 

p= tat +eÂT. (14) 
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Hierin geeft de eerste term rechts de stralingsdruk en ís de 
stralingsconstante a= 7,564 10715 erg em ° KT*, Be tweede term 
rechts geeft de gasdruk en hierin is de gascontante gegeven 
door R= 8,317 107 em° „ Re en is gm het gemiddelde deeltjes= 
gewicht, In een volledig getoniseerd gas met cosmische abundan= 
ties is um # 0,7, Omdat in Röntgendubbelsterren meestal sprake 
is van een geëvolueerd systeem, zal ma groter zijn en stellen we 
verder, pH = 1, zodat we mA niet meer in de vergelijkingen hoeven 
te vermelden en verwarring met het magnetische dipoolmoment is 
uitgesloten, 

Voor de vergelijking voor het energieverlies kiezen we ( 
zie voor de afleiding en andere mogelijkheden Verbunt, 1982): 

á 
Q = EL . (15) 

Hierin is Kk de opaciteit afkomstig van respectievelijk Thomson= 
en vrij=vrij=verstrooiïng: K = 0,4 + 0,66 102 led em°/o. 

Substititie van vergelijking (8) in vergelijking (12) 


geef t: e " 
M= 6nhy. (16) 


Omdat de schijf dun moet zijn, volgt hieruit een ondergrens 
voor de dynamische viscositeitsconstante: n 2 107 g/cm s. Dit 
is vele grootteordes groter den de klassieke waarde als gevolg 
van Coulomb=botsingen (zie bev, Rossi & Olbert, 1970, pe374)e 
De oorzaak van deze grote viscositeit ís echter onbekend, waar= 
door we gedwongen zijn een ad hoc aanname te doen, Zo stellen 
Shakura en Sunyaev, 1973: 


De viscositeit kan bijvoorbeeld veroorzaakt worden door turbu= 
lentie of kleinschalige magnetische velden, In beide gevallen 
moet « £ 1 (zie Verbunt, 1982), Dit 'at-schijfmodel' dient als 


basis voor het model van Ghosh en Lamb (zie G&L I) en de oplos= 
sing în het gebied, dat hierin bekeken wordt (1of £ T$ 10%bom, 
Thomson=verstrooiïng en stralingsdruk zijn hier te verwaarlozen 


), wordt gegeven in Verbunt, 1982: 


hete EN (18a) 
p= 1.3 1072 oct/10 $11/20 „5/8 15/8 jan? (dab) 
Teter 02 er/S #S/10 m/s 4 4 (18c) 
Um BA 1U° „s/° „5/10 8 8, cm/s (18d) 
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hd 





Hierin zijn: ’ pn 
je a en meg, (19) 
Ade 1. 4n 


Een belangrijk facet van de hiervoor behandelde accretie= 
schijfmodellen ís, dat bij de integraties over de z) coördinaat 
zowel Vo als ook V, onafhankelijk van z zijn verondersteld, 
Hoewel deze veronderstelling redelijk lijkt, ís hij zeker niet 
juist: in vergelijking (6) is de afname van de radiële compo= 
nent van de zwaartekracht met toenemende |z| verwaarloosd, De 
correcte vergelijking luidt: 

Eten ì ' 

A ET ke (6) 
Dit geeft dus een verticale snelheidsgradiënt, die een factor 
z/r maal de radiële snelheidsgradiënt ís, Hoewel de verticale 
snelheidsoradiënt dus relatief klein is, geeft hij in het 
schijfmodel van Urpíin, 1983, aanleiding tot merkwaardige effec= 
ten. In het model van Ghosh en Lamb kan de verticale snelheidse 
gradiënt zelfs nog groter zijns we zullen in de volgende para= 
grafen nog op dit probleem terugkomen, Merk echter op, dat ook 
hier de viscositeit een belangrijke rol speelt: als de viscosie= 
teit isotroop is, is dit probleem niet te vermijden; is de vis= 
cositeit niet isotroop, dan kan te veel kleiner blijven dan 
Erp’ zelfs als de verticale snelheidsgradiënt ks wordt, 
Een andere belangrijke veronderstelling in de afleiding 
van de standaard=schijf=vergelijkingen is, dat het radiële 
stralingstransport klein is en de schijf optisch dik is, Veron= 
derstellen we nu, dat de schijf een zwarte straler is, dan 
volgt uit de wet van Stef an-Boltzmann, dat Q evenredig is met 


T Uit vergelijking (13) volgt dan, dat 


„3/4, (20) 


4 
eff ° 
Teff a 
Deze conclusie is volledig onafhankelijk van een eventuele keus 
voor de viscositeit, Waarnemingen van Horne, 1985, geven een 
veel vlakkere relatie voor accretieschijven in cataclysmíische 
variabelen, 

We concluderen, dat nog vele problemen in de standaarde= 
accretieschijftheorie op een goede oplossing wachten, zodat we 
voorzichtig moeten zijn met het doen van quantitatieve uitspra= 


ken op grond van standaard=accretieschijfmodellen, 
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II=2 Het model van Ghosh en Lamb. 


In de inleiding hebben we al een globale beschrijving ge= 
geven van het model van Ghosh en Lamb, We zullen nu aangeven; 
hoe Ghosh en Lamb de structuur van de accretieschijf berekenen 
en hoe hieruit het koppel op de neutronenster kan worden be= 
paald, We moeten nu dus in de schijfvergelijkingen ook rekening 
houden met het magnetische veld en we moeten dus de magneto= 
hydrodynamische vergelijkingen gebruiken, Dat zijn: 

1) behoud van massa: i 


8E - - vlg) (21) 


2) bewegingsvergelijking: 


3) energievergelijking: 
EE + (He) == PV) - Ed - vt A) (25) 
4) de MHD-vergelijking: 
De en 
= Vx(ÙX) — mo 48, (24) 
ef f 


De MHDevergelijking is equivalent met de volgende combinatie 


van de wet van Ohm en de Maxwell-vergelijkingen: 
2 Ira 2 c , 
Î= GO prlE + 7(UXB)) = (VB). | (25) 


Net als bij de standaardeschijf vergelijkingen hebben we 
cok nu weer een toestandsvergelfijkino; een vergelijking voor 
het energieverlies en een vergelijking, die de dynamische vis= 
cositeitsconstante definiëert, nodig, We zullen hiervoor de 
vergelijkingen uit het d=schijfmodel nemen: dat zijn de verge= 
lijkingen (14), (15) en (17) uit de vorige paragraaf, In de 
toestandsvergelijking mag de stralingsdruk verwaarloosd worden, 
zodat met p verder steeds de gasdruk bedoeld wordt, Daarnaast 
hebben we nu ook een vergelijking nodig, die het effectieve gee 
leidingsvermogen, 0 pps definiëert, Net als bij de viscositeit 
is ook voor het effectieve geleidingsvermogen geen simpele keus 
mogelijk, 

Wel kunnen we een eenvoudige parametrisatie voor Ter vin= 
den.(zie G&L I). Daartoe nemen we de radiële component van ver= 
gelijking (25): 5 

J, a Terr (En hi 5 B) nin & 30e). (26) 
De radiële component van het electrische veld volgt uit de ro= 


tatie van de neutronenster en dus van het dipoolveld en is: 
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E add c SsrB,e (27) 
We veronderstellen nu, dat 
B 
dla) = &. | (28) 


Door vergelijkingen (27) en (28) te substitueren in vergelij= … 


king (26) vinden we: 


Ee alef 1 
Terp ankh Fr (Q- Ds) Yo (29) 


Hierin ís de hoeksnelheid (@= va/r) en is Te gedefiniëerd 


door: B B 
ge GD ® Br ao 


Als we nu Ya geven (vergelijk « in de schijfmodellen) îs Toer 
bekend, waarbij we veronderstellen, dat het effectieve geleie= 
dingsvermogen isotroop ís, 

Om nu een oplossing te vinden voor de accretieschijf in 
het model van Ghosh en Lamb te vinden, delen we de oplossing op 
in twee delen (zie ook de inleiding en figuur 5): 

1) de oplossing voor de grenslaag, waarin de magnetische span= 
ning veel groter wordt dan de visceuze spanning, waardoor de 
rotatie snelheid van de materie veranderd van Keplerse rotatie 
naar corotatie met de magnetosfeer, en waarin de materie uit de 
schijf wordt getrokken om via de magnetische veldlijnen naar de 
neutronenster te vallen; de grenslaag ligt tussen Po EP Fo en 
heeft een dikte, 6 « Fo (zie figuur 5). 

2) de oplossing voor het overgangsgebied tussen To en Fr, waar= 
van de structuur bijna gelijk ís aan de structuur van een = 
schijf (zie pagina 12): de visceuze en magnetische spanning 
zijn hier kleiner dan de druk, zodat net als in een standaarde= 
schijf de radiële snelheid klein is en de azimuthale snelheid 


gegeven wordt door Keplerse rotatie, 


De grenslaag, 


Kd 


In de grenslaag mogen we de volgende aannames doen: 
1) axisymmetrie: &5 = 0, 

Dit betekent ook, dat het magnetische dipoolmoment van de 
neutronenster, het impulsmoment van de neutronenster eh het im= 
pulsmoment van de accretieschijf parallel zijn, 

2) de schijf is dun: h& re. 
3) de dikte, Ó, van de grenslaag is klein: $ « re 
Deze eis houdt in, dat er sprake is van een echte grenslaag, 


deWeZe de stroomsnelheid en het magnetische veld veranderen in 
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de radiële richting over een lengteschaal 8» zodat de grens= 
laag een dikte, ó, van de orde van 8, heeft, 

4) de visceuze spanning ís veel kleiner dan de magnetische 
spanninge Omdat B mn neemt de magnetische spanning snel 
toe voor kleinere r,‚ waardoor op zekere afstand de magnetische- 
spanning snel veel groter zal worden dan de visceuze spanning 

en de druk, zodat zich daar een grenslaag vormt, 
5) de schijf is stationair: dr = 0, 
6) het massaverlies uit de schijf (dei, de verticale stroming) 


wordt gegeven door: 


SE M(r) = áge ole)e. (31) 


d 
flux en is g(r) een (onbekende) poortfunctie, die het radiële 


Hierin is Be de geluidssnelheid, M de resterende radiële massa 


prof iel van het massaverlies uit de schijf beschrijft, Ghosh en 
Lamb kiezen hiervoor: 
x) = 1 0 < < 
ol ) 5 > $ X $ Xn (32) 
= (xd Agen) EXE Xe 
Met x,=x,=0,6 en x= (rr ,)/b,- (zie G&L I), 


7) Ye is constant: Ye = Yo’ en van de orde 1. 


Verder wordt dezelfde vergelijking voor het energieverlies als 

in een standaardeschijf (vergelijking (15)) aangenomen en weer 

verondersteld, dat v‚ en Vi niet (of weinig) van z afhangen, We 
kunnen dan de volgende vergelijkingen vinden: 


1) uit het massabehoud: 


4nrhvp = Mi j(r). (33) 
2) de r) component van de bewegingsvergelijking: 
wke. el 1rdp 1 
va == 7*Rr air ee, + ang (9x8) XB). f n (34) 


r 


3) de &) component van de bewegingsvergelijking: 


an TA 2 se 2 . 
d(m GP D/dr = r Cam /dr) + BBF pi (35) 
4) de vergelijking voor hydrostatisch evenwicht: 
2 GM 
p = C ph Es (36) 


N.B. de constante be staat hier, omdat niet aan hydrostatisch 
evenwicht voldaan kan zijn, vanwege de verticale massaf lux, Dus 
moet Es groter zijn dan 3; Er, is nu een vrije parameter van de 


orde 1, 
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5) de energievergelijking (cefe vergelijking (13)): 





- T En J 
Q = kl ze Q = h e (37) 
3keh Terr 
6) uit de d) component van vergelijking (25): 
dB_ _4nG 
z ff 
Tr aes od vB. (38) 
c 


Met de toestandsvergelijking (p = gasdruk) en de vergelijkingen 

(29) en (30) is het stelsel vergelijkingen in principe oplos- 

baar, Ghosh en Lamb lossen de vergelijkingen numeriek op voor 

verschillende waarden van de constanten Yo PP GE, en verschil= 

lende waarden voor de fastnessparameter (zie pag.5)., Als rand= 

voorwaarden kiezen zij: 

1) alle grootheden zijn continu in Fen Fr, en sluiten in Pe 
aan op de oplossing van het overgangsgebied, 

2) in r, is de hoeksnelheid, ® gelijk aan de hoeksnelheid van 
de magnetosfeer, Sys » en wordt B, gegeven door een onge 
stoord dipoolveld: B, = u/r°. 


Enkele conclusies uit hun berekeningen zijn (zie G&L I): 

1) de resultaten zijn vrijwel onafhankelijk van de keuze voor 
de constanten en de fastness, 

2) het B -veld kan door de grenslaag niet volledig gereduceerd 
worden: het resterende veld heeft een grootte van typisch 
0,2 maal het ongereduceerde veld (zie figuur 8). 


3) de positie van de grenslaag wordt gegeven door: 


p= 05 209) oes pf? (zom 7 rd? (59) 


(c.f. pag.2 völ.(1)) 
4) de dikte van de schijf is: Ó = 0,05 T° (40) 





0-4 
a) 
o 4 A ° 1 X 
Figuur 8 De structuur van de grenslaag. voor Ye = 1 en w_=0e3. 
a) de radiële structuur b) de verticale 
structuur 


De dimensieloze grootheden worden in de tekst 
gedefiniëerd, (naar G&L Ii) 
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In figuur 8 is de structuur van de grenslaag geschetst voor een 
typische keuze van de constanten: ke = 24 Je 1; Ws 0,3. De 


hierin gebruikte dimensieloze aenheden zijn als volgt gedefini= 








eerd: a B_(r) 
da SER 
Mm‚(r) 
r- A, ue v/e zh, | 
" ij 8 
4 Een | 
Ö en En % hok Pe En Í 
€ = pen 8 = L. 
Ss s 


De grootheden met benedenindex ‘st zijn achtereenvolgens de 
dichtheid, de druk, de halve dikte.en de temperatuur voor een 


standaard @%=schijf 
Het overgangsgebied, 


De berekening van de structuur van het overgangsgebied 
doen Ghosh en Lamb op dezelfde manier als in de vorige para= 
graaf is gedaan voor een standaardschijf, Eerst worden weer. 
alle aannames gedaan, die we op pagina 10 gegeven hebben, Daare= 
bij geldt alles, wat over de visceuze spanning wordt gezegd, nu 
ook voor de magnetische spanning, Dat geeft ons dan ook dezelf= 
de set vergelijkingen als in de standaardeschijftheorie behalve 
vergelijkingen (7) en (13), waar een magnetische term aan moet 


worden toegevoegd: 





demr? ” 2 ' 
(7) wordt arr + r Te = BBT en (7!) 
4 2 
-_ eat Qt _ 2, dè, 2 J : 
(13) wordt Q = SKER = ybr Gerd) + Toer Fie (13*) 
In vergelijking (7') kunnen we nu de visceuze term weer volgens 


de conventie van Shakura en Sunyaev, 1973, seed juan als t mi =p, Q 
Bovendien moet de magnetische druk kleiner zijn dan de gasdruk, 
omdat de schijf arders uitelkaar wordt gerukt, We kunnen dan de P-| 
totale spanning, Wepes in de vorm Ù er = Op schrijven en we ne= 
men nu aan, dat « constant is over een schaallengte. nr, Verge= 
lijking (7!) wordt dan helemaal gelijk aan de overeenkomstige 
vergelijking uit het @%-schijfmodel,. 
We hebben nu nog een vergelijking nodig, die ons de reduce 
tie van het magnetische dipoolveld geeft, Daartoe schrijven we 


B‚(r) = Hb /r° - B_(r), waarin de eerste term rechts het door 
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het door de grenslaag tot een factor be, (= 0,2) gereduceerde 
dipoolveld van de neutronenster geeft, Deze term geeft in ver= 
gelijking (25) een bijdrage O aan de rotatie-B-term (het is 
immers een dipoolveld). B geeft de verdere reductie van dit 
veld in het OEE B s(Fo ) = OQ en de pine. Tr 
van het overgangsgebied volgt vit 8 s(Fs ) = Hb, 12° . Substitutie 
van deze schrijfwijze voor B, in úeraeld kine (25) geeft: 








dB. 470, kb 
8 eff [e) 
EE SV, 2 (B_ — 5 ). (41) 


We hebben nu alle ingrediënten voor de oplossing behalve 
opnieuw Ye (r). In principe zou Yò gevonden kunnen worden uit de 
®) component van de MHD=vergelijking (vol(24)), maar zolang het 
effectieve geleidingsvermogen onbekend is, is dit niet mogelijk 
en stellen Ghosh en Lamb de volgende bepaling van Yp voor: 

=het azimuthale veld wordt versterkt door de differentiële 
rotatie van de schijf en de magnsetosfeer, dat geef t een 
typische tijdschaal voor de versterking van Ya SSS) » waarin 
Ya een constante van de orde 1 ís, 
„=de dissipatie van dit veld is het gevolg van reconnectie 
met kleinschalige velden in de schijf, wat een karakterise 
tieke tijdschaal Evr/2h geeft, waarin È de dimensieloze recon- 
nectiesnelheid is en v, de Alfvén-snelheid: vj = B9/(4ne)®, 


„vereffening van beide effecten geeft 


- Q 
NS) 


… 
e= P 2h (ane)? (AS zi (42) 
z 


De constante (42/8) volot nu uit de randvoorwaarde, Yol Foor 
Dan wordt gevonden, dat van de orde 0,1 à 1 is, 

Het stelsel vergelijkingen Îs nu dus compleet, De oplos= 
sing van de vergelijkingen voor opnieuw een fastness en 0,3 
is geschetst ín figuur 9, De dimensieloze grootheden zijn weer 
gedefinieerd als in de grenslaag (zie pag. 18), dus ten opzichte 
van de @&-schijf oplossing met dezelfde maar zonder de extra 
dissipatie term ih vergelijking (13'). De dimensiloze straal ís 


gedefiniëerd door y = r/r, en het gereduceerde veld, b door 


_ En 
bout 5 B/ubr 7 
Voor de buitengrens, r‚‚ wordt gevonden, (GEL I val.(41)): 


out? 


r, = 103 7, & -32/73 50/73 Ce, 36/75, (43) 
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5 IO 50 y, 100 


Figuur 9 De structuur van het overgangsgebied.voor w‚=0,5 . 


De halve hoogte €, temperatuur 8, dichtheid Ì en 
het gereduceerde veld bout in dimensieloze eenheden 
(zie de tekst), yer/r, YaPef Fot (naar G&L I) 


Het koppel op de neutronenster, 


Nu we de structuur van de grenslaag en het overgangsgebied 
kennen, kunnen we cok het koppel op de neutronenster berekenen, 
Dit koppel wordt gegeven door de integraal van de impulsmoment= 
flux door een willekeurig gesloten oppervlak, dat de neutronen= 
ster omsluit, Als het systeem axisymmetrisch is,-kan deze inte= 


graal als volgt geschreven worden (zie Lamb, 1977): 


B_e 
K = ba (ret, + Aer u + ned A ds, (44) 


De poloidale snelheid en magnetische veld zijn respectievelijk 


_b 
U_ = V_ë +v 


ä enD =d ë., 
p rr zz p rr 22 


We kiezen nu het oppervlak S, zoals dit is geschetst in 


figuur 10, 








ë 
OUTER 
TEREGION bed TRANSITION —{ 
BOUNOARY ZONE 
LAYER 


_ MAGNETOSPHERIC. 
FLOW 


Figuur 10 Zijaanzicht van de oppervlakken, die gebruikt worden 
om vergelijking (44) uit te werken, 


S eylindrisch oppervlak met straal r en hoogte 2h, 
‘So? 2 vlakken onder en boven de schijf, 
Ss 2 halve bollen in oneindig, 


3 





21. 


De verschillende bijdragen aan het totale koppel worden nu ge= 

geven door (zie GEL II): 

Ko: het koppel van de grenslaag en de materiestroming door de 
grenslaag naar de neutronenster, Omdat op 5, de materiële 
spanning (deie de eerste term onder de integraal in verge=" 
lijking (44)) veel groter is dan de viscueze spanning (cefe 
de derde aanname op pag,10) en de magnetische spanning geen 


component loodrecht op S, heeft, wordt dit koppel gegeven 


4 
doors 
Ki = -pv‚r (r,) 2nr, 2h= 
« 3 (45) 
=M (GMr_)f 2 Ke, 
o o 


Merk op, dat dít koppel gekijk is aan het koppel, dat we in 
de inleiding voor het simpele model vonden (zie pag.3 ver= 
gelijking (2)). 

K‚s het koppel van het overgangsgebied, Er is geen materie= 


stroming door S, en ook de visceuze spanning heeft geen 





2 
component loodrecht op Ss zodat we uit vergelijking (44) 
vinden: 
B_B r 
- zt ee | s 2 2 
Kj-= hs, r( De ) dS = ro alr)B,(r)r dr. (46) 


K_: Omdat Ss in oneindig ligt, verdwijnen alle spanningen op: 53 
en geeft dit oppervlak geen bijdrage aan het totale koppel, 
Als we vergelijking (42) in vergelijking (46) substituören, 


kunnen we K, expliciet berekenen, Sommatie van K, en K‚ geef t 


2 
dan het totale koppel (zie G&L II): 


K = nl. )K (47) 


Hierin is nl) het dimensieloze koppel, dat een functie is van 
de fastnessparameter, Daarnaast is nlw,) afhankelijk van de 
keuze voor de constante Yo in de grenslaag, In figuur 11 ís het 
dimensieloze koppel voor verschillende waarden van Te als funce 
‘tie van de fastness geschetst, We zien ín deze figuur, dat er 
een kritische fastness, Ws bestaat waarvoor het dimensieloze 
koppel nul wordt en van teken verandert, Voor Yorls is WD. 
Analoog aan de afleiding in de inleiding vinden we met 


vergelijking (47) een periodeaf name=lichtkracht relatie: 


(zie Henrichs, 1983) | (48) 
, EN 2, 6/7 
P =(5,9 10" s jr ss Mag MO Re n(w‚) FP Lee 


Als we nu bedenken, dat voor de fastness geldt: 
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DIMENSIONLESS STELLAR ANGULAR VELOCITY 


FiG44.--The dimensionless torque » as a function of the 
dimensionless stellar angular velocity or fastness w‚, for five 
values of the magnetic pitch in the boundary layer, vo. Each 
curve ís labeled with the corresponding value of the pitch. 
Those parts of the curves which continue off the bottom of the 
figure for large w, are shown as dashed curves on the reduced 


le at right. 
scale at rig) (GEL 11) 


f me RE e 
ae Rane 2e et AV en (erS/7) (49) 
Hiermee vinden we dus, dat 
Ld 3 7 8 
B rte 2, (50) 


Kd 
In figuur 12 zijn nu de P's uitgezet tegen de pLS/ re voor een 
aantal Röntgenpulsars, De doorgetrokken líjnen geven de theo= 


retisch voorspelde relaties voor vershillende waarden van het 


log(-P) 


4u 0900-40 





tog(Pi 2/73 


Figuur 12 (G&L II) De periodeafname-lichtkracht relatie, 
De getrokken lijnen geven de relatie als voorspeld 
in vergelijkingen (48) en (50), De labels geven het 


gebruikte dipoolmoment in 1050 G em . 
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magenetische dipoolmoment, zoals we die uit vergelijkingen (48) 
en (49) kunnen afleiden, Merk op dat met name voor Her X-1 (een 
snelle rotator) de overeenkomst veel beter is dan ín de simpele 
theorie (zie figuur 4) het geval was, 

Goede overeenkomst tussen theorie en waarnemingen wordt 
dus voor Her X=1 gevonden als: 
M= 1,3 M_ 30 5 em”, 


3) 
Dat geeft een fastness van ongeveer 0,35, zodat Her X=1 dus een 
fastness heeft, die dicht bíj de kritische ligt, 


Rotatie van de neutronenster met de kritische fastness 


en M= 0,48 10 


geeft een stabiele evenwichtstoestand, De periode die met deze 
evenwichtstoestand overeenkomt, is eenvoudig te vinden door 
substitutie van de kritische fastness in vergelijking (49); dat 
geeft (zie Henrichs, 1983): 


p. & teid es) „el „2/7 mal hea (51) 


e 
Het is dus ook niet uit te sluiten, dat Her X-1 met deze peri= 
ode roteert, 
Deze constatering bracht Henrichs, 1983, tot een alterna= 
tieve verklaring voor de Lange periodeaf name-tijdschaal: 

— als Her X-1 met deze evenwichtsperiode roteert, zal de perie= 
ode toch afnemen als de massa=overdracht en dus L, toeneemt, 
= de toename van de massa=overdracht kan gebeuren op een evolu= 
tionaire tijdschaal (deweze is het gevolg van de evolutie van 

het systeem of de begeleider), dei. Toe 

— uit vergelijking (50) volgt dan onder de aanname, dat de neu= 
tronenster met de evenwichtsperiode kan blijven roteren, een 
tijdschaal voor de periode afname van (zie Henrichs, 1983, en 
Savonije, 1978): 


ee 


T, =Z Tour (52) 


Berekeningen van Savonije, 1978 en 1979, van de evolutie van 
dubbelsterren met Roche-lob=overstroming geven een tijdschaal, 


Tov ge 10° jr, wat-leidt tot een tijdschaal voor de periodeafna= 


me van Her X=1, die goed overeenkomt met de waargenomen waarde, 


T‚=3 10° jr. 





2á 
H-í1 Kritische beschauwing van het model van Ghosh en Lamb, 


In deze paragraaf gaan we verder in op een aantal aannames 
en resultaten uit het model van Ghosh en Lamb, zoals we dat in 
de vorige paragraaf hebben besproken. We zullen hierbij ingaan 
op de door diverse mensen ingebrachte bezwaren tegen het model 
en op de betrouwbaarheid van de quantitatieve resultaten van 


het model, 
Axisymmetriee 


De aanname ván axisymmetrie is gebruikt, om de vergelij= 
kingen niet onnodig gecompliceerd te maken. Een belangrijk ge= 
volg is, dat het magnetische dipoolmoment dan parallel moet 
zijn aan de rotatieas. In dat geval kanzer echter geen sprake 
zijn van een Röntgenpulsar (cef. page2)e 

Toch lijkt het bovenstaande geen groot bezwaar te zijn, 
omdat we niet verwachten, dat de resultaten wezenlijk zullen 
veranderen, als het dipoolmoment een kleine hoek (£ 20°) met 
de rotatieas maakt, Dit introduceert in de vergelijkingen al- 
leen een paar kleine extra termen, die wel verwaarloosd mogen 
worden, Een andere, minder waarschijnlijke verklaring kan zijn, 
dat het magnetische veld vlakbij de neutronenster wordt bepaald 
door hogere multipoolmomenten, waardoor de magnetische pool den 


een aantal graden is verschoven. 
De parametrisatie van het effectieve geleidingsvermogen. 


Zoals de viscositeit de grote onbekende is in de standaard 
schijfmodellen, is het effectieve geleidingsvermogen de grote 
onbekende ín het model van Ghosh en Lamb, Door aan te nemen, 


dat 


1 
Emer orb, en Il=(27) 

B « 
di 11-(29) 


vinden Ghosh en Lamb uit de wet van Ohm en de Maxwell-vergelij= 


kingen (zie pag.14 en 15) het effectieve geleidingsvermogen: 


2 
Cc =| rd _Î 
Jer Fam P Ges) « Il-(29) 


Daarmee is het bepalen van het effectieve geleidingsvermogen 
teruggebracht tot het kiezen van de parameter Yer 

In de grenslaag kiezen Ghash en Lamb voor Ye een constante 
van de orde 1, Yo Omdat in de grenslaag de gasdruk ongeveer 


even groot wordt als B5/en, is deze keuze redelijk; anders 
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wordt de magnetische druk veel groter dan de gasdruk en wordt 

de schijf uiteengetrokken, Het is echter niet uit te sluiten, 

dat Ye in de ans niet constant is, maar bij voorbeeld één 

grootteorde variëert, De invloed hiervan op de grenslaag zelf 

is waarschijnlijk klein, omdat de structuur van de grenslaag 

nauwelijks van de keuze voor Y, afhangt (zie G&L I), De invloed 

hiervan op het overgangsgebied kan echter wel belangrijk zijns 
als Yp in de grenslaag variëert van Yol) = 1 tot YaFo)=10s 
geeft dat een oplossing voor het overgangsgebied, die de= 
zelfde ís, als de normale oplossing bij een grenslaag met 

Yo * 10, In figuur 11 zien we, dat voor zulke grote waarden 

voor Y, het dimensieloze koppel een veel grotere kritische 

fastness geeft dan Yo = 1: dan is w* 0.7 (esfs u * 0,35 

als Á The 

Zolang het werkelijke effectieve geleidingsvermagen in de 
grenslaag onbekend is, geeft dit dus een vrij grote onzekerheid 
in modellen van pulsars die dichtbij hun kritische fastness 
moeten roteren en dus met name in modellen voor Her X=1 , 

In het overgangsgebied wordt Yò berekend door de verster= 
kingstijdschaal van het azimuthale veld (bepaald door de diffe= 
rentiële rotatie van de schijf en de magnetosfeer) en zijn dif= 
fusietijdschaal (bepaald door reconnectie) ek te stellen. 
Dat geeft (zie pag.19): 


ke, 


- de 2h Cano)? (== II=(42) 


waarin oi gevonden kan worden uit de voorwaarde, dat Y 
continu moet zijn in Tr, deweZe Yol Fo ) = Yo* Wanneer we dat 
doen, vinden we sden voor ë van 0,1 tot 1 ele is van de orde 
13 zie GEL I), Dat betekent, dat de reconnectie moet gebeuren 
met snelheden in de orde van de Alfvén=snelheid, Het is echter 
zeer de vraag of dergelijke hoge reconnectiesnelheden wel moge= 
lijk zijn: de geometrie van dit probleem wijkt wel sterk af van 
de gebruikelijke geometrie in X=type reconnectie en dat zou de 
reconnectie snelheid sterk kunnen reduceren (zie bev. “::……. 
Vasyliunas, 1975), 

Bavendien ís hat zeer de vraag, of dergelijke reconnectie 
processen wel beschreven kunnen worden met een effectief gelei= 
dingsvermogen (zie Vasyliunas, 1979), Daarnaast wordt bij X- 
type reconnectie de magnetische energie omgezet in beweging van 


de materie, terwijl in de Keplers roterende schijf de beweging 
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beperkt ís, 
De stroming uit de schijf, 


In de grenslaag wordt de materie in de schijf gedwongen om 
langs de veldlijnen te gaan bewegen en dus de schijf invertica= 
le riehting te verlaten, Hoe dit gebeurt ís in het model van 
Ghosh en Lamb onbekend, omdat de structuur van het magnetische 
veld boven de schijf nog onbekend is, Daarom modelleren zij het 
massa verlies met de poortfunctie g(r) (zie pag,16), Zij vinden 
dat de structuur van de grenslaag ongevoelig is voor de precie= 
ze keuze van g(r). 

Een belangrijk probleem wordt echter aangegeven door Anzer 
en Börner, 1983, Zij menen, dat als gevolg van de reductie van 


het magnetische veld met de factor sân de grenslaag het resul= 


5 
terende veld sterk afwijkt van een dipoolveld: in figuur 13 is 
het totale veld geschetst voor het geval dat 2 deel van de flux 


door de opening in de schijf (r = ro) gaat en slechts ï deel 


door de schijf gaat, We zien hier, dat langs de veldlijnen, die 





Fig! 4. Aligned rotator with t of the flux through the disk 
(Neutron star not drawn to scale) 


door de grenslaag gaan, de potentiële zwaartekrachtsenergie een 
minimum heeft, De materie die langs de veldlijnen naar de polen 
van de neutronenster stroomt moet dus over een forse potentiaal 
berg, wat haast onmoaeli jk lijkt, 

Zoals we hierboven al opmerkten, is de werkelijke velde 
structuur onder en boven de schijf eigenlijk onbekend (zie ook 
G&L I). Het is dan ook niet uitgesloten, dat het veld in zijn 
oeheel door de accretiestroming naar de delen wordt samenge= 
drukt, zodat het resulterende veld zowel binnen als buiten Fo 


dipool=achtig blijft met een dipoolmoment van respectievelijk 
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kp en Ân en de potentiële zwaartekrachtsenergie in het equator= 


vlak steeds een maximum heeft, 
De resultaten voor Her Xe 


In zowel de verklaring van Ghosh en Lamb (zie G&L II) als. 





die van Henrichs, 1983, voor de langzame periodeafname van Her 
X=1 wordt verondersteld, dat Her X-1 roteert met een snelheid, 
die bijna gelijk is aan de kritieke fastness, Met de standaard 


keuze voor de grenslaagparameters betekent dit, dat we, “ 0,35, 





Dat geeft tenslotte een magnetisch dipoolmoment voor Her Xe1 
van ongeveer 0,5 10% 5 em° (zie ook pag.23). 

In het spectrum van de Röntgenpulsen van Her X=1 wordt een 
lijnstructuur waargenomen (zie Trümper et al, 1978, en figuur 


13). Trümper et al, 1978, verklaren deze structuur als een ge= 


N 
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Figuur 13 Het Röntgenspectrum van de Her X=1 pulsene 
De getrokken lijn geeft de beste aanpassing met een expo= 
nentieelspeectrum en een Gaussische emissielijn bij 58 keV, 
De onderbroken lijn geeft het totale Röntgenspectrum van 
Her X=1 waargenomen door 050=8 ín 1975, Merk de goede over 
eenkomst tot 40 keV op, (uit Trümper et al, 1978) 





28 


quantiseerde electronen=cyclotron=emissielijn bij 58 keV, 
De gequantiseerde energieniveaus voor electronen in een 
homogeen, sterk magnetisch veld worden gegeven door: 


Pp 
Es = me Kk + (2)? + (2j +s + ne} = 1}, (1) 


waarin j = 0,1,2,eee en $ = u respectievelijk het impulsmo= 
ment= en spin=quantumgetal zijn, B de magnetische veldsterkte 
de, Bm (mo fak) =4ä 4072 


Cr 
impuls langs het magnetische veld is, Trümper et al, 1978, ver= 


G, m de electronmassa îs en p‚ de 


klaren de emissielijn bij 58 keV als afkomstig van de overgang 
van het eerste aangeslagen niveau (j= 1, 5 = =1 of j = 0, e=1) 
naar de grondtoestand (j = 0, s = =1)e Als we nu p‚ in verge= 
lijking (1) verwaarlozen, vinden we een magnetisch veld in het 
gebied, waar de lijn gevormd wordt van 5,3 1012 G, De tweede 
harmonische zou dan bij 110 keV moeten optreden en in figuur 13 
zien we, dat bij 110 keV inderdaad weer een piek te zien is. 

Een andere mogelijke verklaring voor de líijnstructuren in 
dit spectrum is absorptie met de eerste absorptielijn bij 40 
keV en een tweede bij 80 keV, Net als hierboven vinden we uit 
vergelijking (1) het magnetische veld in het gebied waar, deze 
lijnen gevormd worden, dei, B= 4 1012 G. Op theoretische grone 
den verdient de emissielijn=verklaring echter de voorkeur (zie 
Basco & Sunyaev, 1975), 

Algemeen wordt aangenomen, ‘dat de harde Röntgenstraling 
afkomstig is van het oppervlak van de neutronenster, Er moet 
dan gecorrigeerd worden voor de gravitationele roodverschuiving 
en Trümper et al, 1978, geven dan een magnetisch veld op de 
polen van Her X=1 van 5,8 tot 7,3 1012 G, Dat betekent, dat het 
magnetische dipoolmoment van Her Xe1 ongeveer 3 10°? G em“ ise 

Het magnetische dipoolmoment, dat we zo vinden, is dus veel 
groter dan het magnetische dipoolmoment, dat volgt uit het moe 
del van Ghosh en Lamb, Een mogelijke (alhoewel niet erg waar= 
schijnlijke) verklaring is, dat dichtbij het oppervlak van Her 
Xt hogere muitipaolmementen belangrijk worden, 

Ook kunnen we voor Her X=1 een relatie tussen de fastness 
en het magnetische dipoolmoment afleiden uit de vergelijkingen 
Il-(39) en Il=(49), Met P= 1,24 5, M= 1,4 MQ en m= 1017 o/s 


(zie bev. Henrichs, 1983) vinden we: 


Pe 1e7 10°0 / Gems (2) 
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Daarnaast zien we in figuur 11, dat voor grotere waarden voor 
de constante Yo ook de kritische fastness behoorlijk groter 
wordt, Voor te 8 bijvoorbeeld zien we, dat de kritische faste 
ness ongeveer 0.7 wordt, Uit vergelijking (2) volgt den, dat 
het magnetische dipoolmoment ongeveer 1 105 G em” is; een Te=: 
sultaat dat al iets beter is dan het resultaat van Ghosh en 
Lamb. We zien echter ook uit vergelijking (2), dat in dít model 
het magnetische dipoolmoment nooit groter mag zijn dan 1,7 10° 
G em, Voor grotere dipoolmomenten wordt de fastness groter dan 
1, wat betekent dat er geen accretie meer plaats vindt, omdat 
op de grens met de magnetosfeer de materie zal worden wegge= 
slingerd, 

Er blijft dus een belangrijke discrepantie tussen de waars= 


nemingen en het model van Ghosh en Lamb bestaan, 
Het verticale, hydrostatische evenwicht, 


In het model van Ghosh en Lamb wordt de druk in het equa= 
torvlak gegeven door dezelfde vergelijking als in de standeard= 


schijfmodellen (zie pag.11): 
p= den” 5, Ir-(9) 
r 


Boven en onder de schijf moet de druk kleiner zijn, zodat de 
drukgradiënt de kleine verticale component van het zwaarte= 
krachtsveld kan compenseren (zie pag.11). 
De magnetische druk ín Ben En wordt gegeven door: 
B 

Er 0de dt 
Ook deze moet dus kleiner zijn dan de druk als gegeven in ver= 
gelijking Il-(9), In grote delen van het overgangsgebied zal 
be veel groter zijn dan 1, zodat we uit het bovenstaande de 


volgende ongelijkheid kunnen afleiden: 
202 2 GM 2 « 
Yarr Sten” G= ben (4) 
® z se 


In deze vergelijking substituëren we nu de vergelijking II=(42) 


voor Se in het overgangsgebied met (fa/8) a 13 we vinden dan: 
2 2 22 
eh (ee, = Bs) < Zeh %. (5) 
Buiten de corotatiestraal zal SN snel veel kleiner worden dan 


is 


daan. Met andere woorden: de oplossing is instabielë 


en kan daar niet aan de bovenstaande ongelijkheid zijn vol= 
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In feite is het bovenstaande probleem een gevolg van de 
slechte behandeling van de verticale structuur in de standaard= 
schijfmodellen, Daarom zullen we nu de relevante vergelijkingen 


opnieuw afleiden: 
De bewegingsvergelijking. 


We schrijven de bewegingsvergelijking in het met de magne= 
tosfeer coroterende stelsel, Dat geeft in het stationaire geval 


de volgende vergelijking: 


HIE + Bd + Sexo) = 
B? pe Er an (6) 
— pVÍ - Vr tam = Wet + zr. 
Met de gebruikelijke aannames uit de standaard=schijftheorie 
(zie pag.10) en door alleen de termen met het radiële magneti= 
sche veld te verwaarlozen vinden we nu: 


r) component, 2 B 


Vv 
mt: 2 EN? 
Et Bisto + MST = 2 + ampt + 


In het equatoervlak zal de laatste term te verwaarlozen zijn, 


2 
ek. “lr 


es 
ie) 


maar daaronder of daarboven kan eg zo klein worden, dat deze 
term wel belangrijk wordt, We vinden dus, dat de rotatie in het 
equatorvlak Keplers is, maar dat de rotatie onder en boven dit 
vlak belangrijk van Keplerse rotatie kan afwijken, Dit is ook, 
wat we verwachtten, omdat ver onder en boven de schijf de mate= 
rie moet coroteren in de magnetosfeer, Deze laatste term ver= 
telt ons dan ook waar de grens tussen de schijf, waarin de ma= 
terie Keplers roteert, en de magnetosfeer, een starre rotator, 
ligt, 

®) component, 


B 

zev‚Sr(rvg) ze GN, + „ze are (a) 
Als de afwijking van de Keplerse rotatie, zoals die uit de vo- 
Tige vergelijking volgde, nu maar klein ís zolang gp nog onge- 
veer gelijk is aan de dichtheid in het equatorvlak, mogen we 
deze vergelijking op de gebruikelijke manier over z integreren 
en vinden we dezelfde vergelijking als in het model van Ghosh 
en Lamb, deie vergelijking Il=(7!), 


z) component. 


dz tE EE Ca} 
T 
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Dit is nu de vergelijking voor hydrostatisch evenwicht, Hierin 
is p_ de gasdruk, ín het overgangsgebied is de stralingsdruk te 


verwaarlozen (zie G&L I)e 
De magnetische vergelijkingen, 


In het met de magnetosfeer coroterende stelsel mag het e= 





lectrische veld in de wet van Ohm gelijk aan nul verondersteld 
worden; de MHD=evergelijking verandert niet, Deze laatste geeft 


voor de r en $ componenten van het magneet veld de volgende 





vergelijking (zie pag,14): 


2 2 | 
B nn EN Je (10) | 


(v - vB) = (B 

92 zZ T‚& T‚Ò z ANO er dz2 Tr, 

Omdat voor de Tr en & componenten dezelfde vergelijkingen gelden, 

mogen we de volgende veronderstelling doen: 
v_(z) v‚(z=0) B‚(z) 
Er = = i (11) 
va lz) vg (2-0) BCZ) 

De vergelijking voor B, kan ongewijziad uit het model ven 


Ghosh en Lamb worden overgenomen, deie vergelijking Inl 41), 
De energie=vergelijkinge 


Ook de energie vergelijkingen nemen we onveranderd uit het 


model van Ghosh en Lamb over, dei, vergelijking II=(13'), 


Met dezelfde keuzes voor de toestandsvergelijking, de vis= 
cositeit,. het energieverlies en, het effectieve geleidingsvermoe= 
gen als in het model van Ghosh'en Lamb is het stelsel vergelij= 
kingen weer compleet, | 
We construêren nu een nieuwe oplossing voor de verticale 


structuur van de schijf, Daarbij doen we de volgende aannames: 


1) in het equatorvlak wordt de oplossing gegeven door de stan= 

daard=schijf oplossing met of= 1, Deze is gegeven in verge= 
lijkingen II=(18b t/m d). Deze aanname kunnen we rechtvaardigen, 
als het grootste deel van de materie Keplers roteert (delle Zie 
als de grenslaag tussen de schijf en de magnetosfeer onder en 
boven de schijf dun is of als de dichtheid in die laag veel 
kleiner is dan de dichtheid in het equatorvlak), zodat we voor 
vergelijkingen (7) en (8) de overeenkomstige vergelijkingen uit 
het model van Ghosh en Lamb mogen gebruiken, en omdat de struc= 
tuur van het overgangsgebied in het model van Ghosh en Lamb 
niet sterk afwijkt van de standaard=schijfoplossing (zie figuur 
9), De afwijking blijft steeds ruim binnen één orotteorde, Dit 
geldt ook voor De en we zullen B, = 0.2 m/e° stellen, 
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2) omdat we hier niet geïnteresseerd zijn in het verticale e= 
nergietransport, maar willen weten hoe het magnetische veld 

de schijf globaal verandert, nemen we ean dat de schijf in ver= 

ticale richting isotherm ís met de temperatuur als gegeven in 


vergelijking Il=(18c), 
3) de toestandsvergelijking wordt gegeven door: 
pj RT (het gem, deeltjesgew, is 1; zie paa.,12) 
4) het effectieve geleidingsvermogen is hetzelfde als in het 
model van Ghosh en Lamb, Dit vinden we door substitutie van 
vergelijking II=(42) in vergelijking II=(29): 


G 


ef 5 e Carp). stg « (12) 


z 


5) in het overgangsgebied mogen we buiten de corotatiestraal 
in de laatste term het deel met B° verwaarlozen, omdat B? 
hier maar zwak afhankelijk is van r (uit vergelijking II=(42) 


volgt met substitutie van vergelijking II-(18a en b) BÔmr). 


Daarmee is nu op elke plaats, r‚ de volgende set vergelij= 


kingen oplosbaar: 





vz B2 

ó 2__ GM ò 

rt 25 * Ast = 2 *Gner ’ (7) 
2 

] En en 

Sz(p, * En zes | en (2) 


2 


ò 9 
Sr(voBp = 42) = GPR ie (10) 


N.B, we verwaarlozen B en dan wordt de z afhankelijkheid van 
Be dus helemaal bepaald door Bs omdat B bij benadering onaf= 
hankelijk is van zeals z klein is, 


We kunnen nu vergelijkingen (9) en (10) herschrijven als 


integraal vergelijkingen: 
z B 
A GM p) 
e(z) = RT Á (zes + zE 57E ez + e(z=0). (13) 
o T 
F4 SNT pp Ô 
B (zZz) = —_l VB B_)dz. 14 
ple , Sung - venen (14) 
In deze laatste vergelijking is de volgende randvoorwaarde 
gebruikt: 
ANT op 








A38 


GOE 





Ne 


004 


GG) 


De vergelijkingen (7), (11) en (12) met als randvoorwaarden de 
oplossing voor de @-schijf in z = O0 ‘kunnen nu numeriek worden 
opgelost, We gebruiken hiervoor het IBM standaard programma 
CSMP III (Continuous System Modelling Program), dat geschikt is 
voor het oplossen van een stelsel niet-lineaire, gewone diffe=. 
rentiaalvergelijkingen, Deze berekeningen werden gedaan op de 
DEC=20 van het Computer Instituut Woudenstein te Rotterdam, 

In appendix A geven we een korte beschrijving van CSMP III 
en in âppendt% B beschrijven we het specifieke CSMP-programma, 
waarmee we het bovenstaande stelsel vergelijkingen hebben opge= 
lost. 

In fiouur 14 geven we een karakteristiek voorbeeld voor de 
bepaalde oplossingen, We zien hier, dat vg inderdaad alleen in 


een dunne grenslaag afwijkt van de Keplerse rotatie, 
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z/h 
Figuur 14 De verticale structuur voor r = ST en 0 = 0,30, 
Gamma is hier gedefiniëerd door: B‚(2)/8,E Y(z). 


De hier gebruikte grootheid h ís niet de halve hoogte van de 
schijf, zoals deze gedefiniëerd is in de standaard=schijftheo= 
Tie, maar is de schaalhoogte van de druk, zaals die volgt uit 


vergelijking (9) met ip 0 
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(16) 
Deze verschilt slechts een factora/2 van de h uit de standaarde 
schijftheorie, | 

We zien ook in figuur 14, dat de grens met de magnetosfeer 
nog beneden z = h ligt. Als deze grens te laag wordt is echter 
met massatransport in de schijf. niet meer constant, Om er nu 
voor te zorgen ‚ dat het massatransport in de schijf wel overal 
gelijk is, vermenigvuldigen we de dichtheid in het equatorvlak 


met de factor C, zódat Cz = h (z is de grens tussen de 
max max 


a 
schijf en de magnetosfeer). 
Uit de zo verkregen resultaten lezen we dan de waarden van 


year in z dat geeft ons dan de Yb in ons model, De zo ge= 


max? 
vonden waarden voor Kn gegeven ín figuur 15 voor de typie= 


M- uu Mo 
id M pd oi q/s 


W: 0.30 





gecorrigeerd voor 
p 3 druk - evenwicht 






1 r 5 10 —_p P/N 5o 400 


Figuur 15 De resultaten voor vg rekeninghoudend get verticaal 
drukevenwicht: karakteritieke waarden voor magsa en 
accretiesnelheid en w‚= 0,30, 


de onderbroken lijn geeft de oplossing voor Yò in vergelij= 
king IIl=(42) volgens G&L I-.en II, 


sche neutronenstermassa, M = 1,4 Mos „=accretiesnelheid, M=1017 


g/s, en fastness, W, = 0,30, Duidelijk is, dat de door ons ge= 
vonden waarden voor Ye veel kleiner worden dan de waarden, die 


Ghosh en Lamb geven voor Vet 
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Een consequentie van deze veel lagere waarden voor Yb is, 
dat ook de bijdrage aan het totale koppel van het gebied buiten 
de corotatiestraal veel kleiner wordt, Het is juist dit gebied, 
dat een negatieve bijdrage aan het totale koppel geeft (zie pag. 
6 en 7). Het verschil wordt bovendien steeds groter als Tr groe 
ter wordt en volgens Ghosh en Lamb is de bijdrage aan het nega= 
tieve deel daar juist het grootste 

We concluderen hieruit, dat de fastness, waarvoor het kop= 
pel op de neutronenster nul wordt, veel groter is dan Ghosh en 
Lamb voorspellen,of dat het koppel voor geen enkele fastness, 
kleiner dan 1, nul wordt,-maar steeds positief blijft. 

Daarnaast werd de constante C in een aantal berekeningen 
ongeveer 50 & 100, zodat de schijf daar wel erg veel dunner en 
dichter zou moeten zijn dan een standaardeschijf, Het is de 


vraag of dit wel reëel is, 


Uit al deze problemen wordt duidelijk, dat op grond van 
het model ‘van Ghosh en Lamb eigenlijk alleen maar de qualita= 
tieve uitspraak gedaan kan worden, dat er wel eens een critie 
sche fastness zou. kunnen bestaan waarvoor het totale koppel nul 
wordt, Voor verdere quantitatieve uitspraken zijn betere bere. 
keningen van de verticale structuur van de schijf en betere be= 
palingen van het effectieve geleidingsvermogen nodige 

Het is bijna onmogelijk om op theoretische gronden het ef= 
fectieve geleidingsvermagen te. bepalen; ook in de laboratoriume 
plasmafysica ís men meestal niet in staat op zuiver theoretie 
sche gronden de waargenomen transportcoëfficiënten te verklaren, 
Zowel de viscositeit als het effectieve geleidingsvermogen moe= 
ten în deze schijf modellen vele ordes van grootte van hun 
klassieke waarde, die het gevolg is van de Coulomb=botsingen, 
verschillen, We verwachten dan ook, dat beide het gevolg zijn 


van dezelfde turbulente of andere verschijnselen in de schijf. 


« 
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MH-2 Alternatieven, 


In de vorige paragraaf hebben we laten zien, dat de pro= 
blemen in het model van Ghosh en Lamb veroorzaakt worden door 
de veronderstellingen omtrent de verticale structuur en het ef= 
fectieve geleidingsvermogen van de accretieschijf, Als we het 
model van Ghosh en Lamb willen verbeteren moeten we ons dus op 


deze twee punten richten, 
De verticale structuur. 


Naast het al uitvoerig in de vorige paragraaf besproken € 
hydrostatische evenwicht is een andere belanrijke veronderstel= 


ling met betrekking tot de verticale structuur, dat 


98 8 
& & 
er Takeh d II-(28) 


We zien hiervoor echter geen directe aanleiding en we kunnen 
met evenveel recht veronderstellen, dat B verandert in een dun 
laagje, met een dikte d veel kleiner dan de halve dikte van de 
schijf (d & h), aan de onder= en bovenrand van de accretie- 
schijf (zie figuur 16). Dan is duss 

| ed (17) 

z d 

Merk op, dat er nu alleen in dit dunne laagje een radiële 
stroom loopt en dat in de rest van de schijf de radiële stroom 
verdwijnt, omdat (vB), = 0 buiten deze dunne laag. De radiële 
component van de wet van Ohm wordt dan (in de rest van de 
schijf): 
Em ie Br. (18) 





z=h 
Î A 
| 
Ì 
hy le | 
! 
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z=0 4 


Figuur 16 De verticale structuur van het magnetische veld en 
de stroming in de schijf. 
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Dat houdt ook in, dat de veronderstelling in vergelijking I= 
(27) niet meer juist ís, 

Verder worden de verticale gradiënten in het laagje tussen 
de schijf en de magnetosfeer heel groot, We veronderstellen nu, 
dat in dit laagje alle andere termen in de vergelijkingen ver= 
waarloosd mogen worden ten opzichte van de verticale afgelei= 
dene 


De & component van de bewegingsvergelijking wordt dans 


2 
0 =D) - EO (19) 
De fg component van de MHD-vergelijking wordt dan: 
2 2 
d UNE NE kf 
zt pez) TT AGopr ye): (20) 


We gebruiken hierin weer het met de magnetosfeer coroterende 
coördinaten systeem, De azimuthale snelheid moet in het laagje 


dus nul worden, zodat ook 


Òv v 
erder s an 


Als we nu vergelijkingen (17) en (21) ín de vergelijkingen (19) 


en (20) substituëren, vinden we 


2 

de Sn en (22) 
B Orr 

B = v (efo Ìs (23) 
$\° Terr | 


Om onze veronderstelling, dat d veel kleiner ís dan h, te kon= 
troleren voeren we de het Hartmann=getal, M‚, in (cef… Jackson, 
1962, page475=479) 


2 
8_T „‚h 
n= l- (HEt, (24) 


d e2n 
Als de veronderstellingen en de resultaten consistent zijn, 
moeten we dus een Hartmann=getal vinden dat veel groter is dan 
1e 

Een probleem is echter, dat we het effectieve geleidings= 
vermogen en de dynamische viscositeitsconstante niet kennen. 
Voor de laatste kiezen we de waarde, zoals die volgt uit de verm 
Gelijkingen Il-(8) en Il-(17) voor een schijf (zie Verbunt, 
1982): | 


np 


pese. | Oe (25) 


6) 
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Ter illustratie kiezen we voor het effectieve geleidingsvermo= 
gen twee verschillende mogelijkheden: 

1) het klassieke geleidingsvermogen (zie Rossi & Olbert, 1970, 
pag.374): 


D= 7 107 r/ una)! dl, (26) 


Hierin is LnA de Coulomb-logaritme, deze Ís van de orde 1 à 10, 
2) een turbulent geleidingsvermogen (grootste cellen in turbu= 


lentie zijn ongeveer 0,1 h)s 


T & 10° se ( Nak (27) 
Et un ‘Cs . 


Hierin is c, de geluidssnelheid, ce, =(p/e). 


In de bovestaande vergelijkingen substituêren we nu de 
resultaten voor het «-schijfmodel van vergelijking II-(18) en 
we vinden de volgende resultaten: 
bij klassiek geleidingsvermogen 


2 4010 “3/10 41/20 „1/8 9/8 d 


Mm 


4 


E z’ 


32 go1/10 „13/20 3/8 ‚1/8 


B m 


® 


en bij turbulentie 
15 q1/20 4=17/40 m11/16 „21/16 je 


R 


4 10 (R — Rs) 


Mm & 2 10 Eis 
B (A5) 


we zien hierin, dat inderdaad het Hartmann=getal kleiner is dan 


IR 


da ge/20 211/40 „1/16 „1/16 


1, maar de magnetische velden worden zó groot, det de magneti= 
sche druk weer vele orden groter wordt dan de gasdruk in de 
schijf (cef. pag,29), Ook deze methode brengt ons dus niet 


dichter bij een oplossinge 
Het effectieve geleidingsvermogene 


Zoals we al in de vorige paragraaf opmerkten, heeft het 
effectieve geleidingsvermogen waarschijnlijk dezelfde oorsprong 
als de viscositeit, We hebben dus een schijfmodel nodig, dat 
een expliciet model geeft voor de oorsprong van de viscositeit 
en waaruit bovendien een waarde voor het effectieve geleidingse 
vermogen kan worden afgeleid, Voor zover ons bekend is, is het 
enige schijfmodel dat aan deze voorwaarden voldoet, het model 
van Ichimaru, 1977, 

In de volgende paragrafen zullen we de belangrijkste ele= 


menten van dit model kort weergeven om daarna de berekeningen 
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van Ghosh en Lamb te herhalen met als basis dit schijfmodel van 


Ichimaru în plaats van het of-schijfmodel, 


Daarna gaan we ín op andere alternatieven, waarin het pro= 
bleem van de langzame periodeaf name van Her X=1 op volledig an= 
dere wijzen wordt verklaart, Dit betekent, dat met name de ver= 
onderstelling, dat het magnetische dipoolmoment parallel is aan 
de rotatieas, wordt verlaten. In het model van Anzer en Börner, 
1980 en 1983, bijvoorbeeld staat het magnetische dipoolmoment 


loodrecht op de rotatieas, 
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IVe1 Het schijfmodel van Ichimarue 


In het accretieschijfmodel van Ichimaru, 1977, is het de 
magnetische turbulentie, die door de differentiële rotatie van 
de accretieschijf wordt opgewekt, die voor het impulsmoment= 
transport zorgt en daarbij ook een anomaal effectief gelei= 
dingsvermogen geeft, In Ichimaru, 1975 en 1976, worden het ef= 
fectieve geleidingsvermogen bij de aanwezigheid van magnetische 
turbulentie, de r@=component van de spanningstensor, die het 
gevolg is van de magnetische turbulentie, de magnetische ener= 
giedichtheid in de turbulentie en het energiespectrum van de 
magnetische turbulentie afgeleid voor het geval van een diffe= 
rentiëel roterende schijf, 

Als de turbulentie isotroop ís, wordt het volgende effec= 


tieve geleidingsvermogen gevonden (Ichimaru, 1976, vgl.(16)): 


+ 2 
_ Snm{T 1E Es 71 
Spe © en : CE % Ee 9. (1) 


Hierin is n de deeltjesdichtheid van de ionen, T de temperatuur 
in erg's (deweze T = T(K)k, ), m de massa van de ionen en is 

(B kt" B> de gemiddelde mes heldeslk energiedichtheid in de tur= 
je velden met golfgeteal k, Om dat de electronen sterk aan 
de magnetische velden gebonden zijn is in vergelijking (1) al= 
leen de bijdrage van de ionen aan de stroom meegenomen (zie 
Ichimaru, 1975). e 

Het energiespectrum van de magnetische turbulentie wordt 


gegeven door (zie Ichimaru, 1976, val.(18)): 


3 2 
_ 2n,” 68°) 
B erm) (2) 
el 4 Ge Id 
waarin <58) = 3 <B_-B> de gemiddelde totale (3) 


energiedichtheid in de magnetische turbulentie is en waarin L 
de maximale golflengte in het turbulentiespectrum is. 

De totale energiedichtheid in de magnetische turbulentie 
wordt voor een differentiëel roterende schijf gegeven aen (zie 


Ichimaru, 1976, val.(27)): 


Ee ie le Le, ea (4) 


Bn nT 


2 
waarin ve zn? at (zie Ichimafu, 1976, vol.(22}).(5) 
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We substituëren nu de vergelijkingen (2) t/m (5) in verge= 
lijking (1) en vinden dan het effectieve geleidingsvermogen: 
—_Î 


3 -1 
3 2 (Lt pT. 1 an 
Sefr Eik eG) le È (25 in 


It 


(6) 


4 


2 1 
Zok 2 21 AR 
Die lr ’ 
Hierin hebben: we gebruikt, dat de som in deze vergelijking on= 


geveer gelijk is aan de grootste term in de som: 
L 


1 4 
K k° k 5 2n 
min 


Op soortgelijke wijze kunnen we de gemiddelde rg=component 
van de spanningstensor vinden (Ichimaru, 1976, vol,(15)): 
1 T 
hs cher Ks ë (7) 
In een differentiëel roterende schijf kan nu afgeleid worden, 
dat (zie Ichimaru, 1976, vgle Ee 


ee (/mt » &, (8) 


Als we deze vergelijking in vergelijking II-(8) substituëren, 
vinden we de volgende definitie voor de turbulente viscositeit: 
Td 
Sr In (9) 


Samen met vergelijking (8) vinden we nu: 


n= en daj ks (10) 
Er) 

Met deze viscositeit en de standaard=schijf vergelijkingen 
uit hoofdstuk II kan een accretieschijfmodel gemaakt worden. 
In deze modellen is de grootste golflengte, Lt, in de magneti= 
sche turbulentie een vrije parameter, vergelijkbdar met de « 
uit de «=schijfmodellen, Als we voor L de halve dikte van de 
schijf, h, kiezen, vinden we het schijfmodel van Ichimaru, 1977, 
Bedenk wel, dat in de enn ae nu ook de magneti= 
sche druk ín de turbulentie, deie <ÁB 2/an, moet worden meege= 
nomen. 

In het voor ons interessante gebied rond r & 10° em zijn 


de stralingsdruk en de vrij=vrij=verstrooiïng te verwaarlozen 
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en wordt de oplossing voor de standaard=schijfvergelijkingen 


gegeven door (zie Ichimaru, 1977, vol.(47) en (48)): 


h = 6,0 1075 a1/ „7/20 „11/20 cm (11a) 

= 1,4 1010 #2/5 „11/20 „33/20 fam” (11b) 
T Sat doe a2/° „5/10 „5/10 © (11e) 
ús dez 10° a2/° IR „2/5 cm/s. (11d) 


T 
‚Met dit schijfmodel kunnen we het model van Ghosh en Lamb 
opnieuw doorrekenen, We bekijken eerst de grenslaag en verge= 
lijken het effectieve geleidingsvermogen uit vergelijking (EN 
met het effectieve geleidingsvermogen in de grenslaag van het 


“model van Ghosh en Lamb (zie G&L I), deie vergelijking II-(29): 


2 2 
Ee“ ele ef _ 2 (2 kor! 
Ba mT ot Bnr CP) ER ee MO 


In het model van Ghosh en Lamb is Yo Van de orde 1 en de dikte, 
8, van de grenslaag is ongeveer van dezelfde grootteorde als he 
Hieruit mogen we concluderen, dat beide modellen ongeveer hete 
zelfde effectieve geleidingsvermogen geven in de grenslaag. 
Omdat de structuur van de grenslegg ín het model van Ghosh en 
Lamb nauwelijks van de voor Yo gekozen waarde en dus nauwelijks 
van de precieze waarde van het effectieve geleidingsvermogen 
afhangt, mogen we veronderstellen, dat we in het model met het 
effectieve geleidingsvermogen uit Ichimaru, 1976, een grenslaag 
zullen vinden met dezelfde karakteristieken als een grenslaag 
in het model van Ghosh en Lamb (zie G&L I en page17). 

We veronderstellen dus, dat de grenslaag dun is, é «rs 
dat de grenslaag het magnetische dipoolveld van de neutronen= 
ster slechts met een factor 0.2 kan reduceren en dat de straal, 


waar de grenslegg begint, gegeven wordt door 
r,= 0,5 „(0) = 0,5 tr (zem)=1/? m2/7, II-(39) 


De structuur ‘van het overgangsgebied berekenen we op de= 
zelfde manier, zoals die op de pagina's 30 t/m 32 uiteen is ge= 
zet voor de berekeningen in het model van Ghosh en Lamb, die 
gecorrigeerd zijn voor het verticale drukevenwicht, Dat houdt 
in, dat we veronderstellen dat de radiële structuur van de 
schijf niet of maar weinig verandert door de interactie met het 


magnetische veld van de neutronenster (c.f‚… page31). 
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Wel berekenen we opnieuw de radiële structuur van het mag= 
netische veld, die ons onder andere de uitgestrektheid van het 
overgangsgebied geeft, Hiervoor gebruiken we vergelijking II- 


(41), die ons het verticale veld geeft (zie pag.19): 


dB ang, bu 
s eff P fe) 
 Ú En (B; = TE ). Il(41) 


We substituëren hierin vergelijking (11d), met m= @ = 1, en 
vergelijking (6), maar nu met L = 4h, Dit laatste wordt gemotie 
veerd door het feit, dat het geleidingsvermogen sterk van L af= 
hangt: Als we L kleiner kiezen wordt B veel groter en we ver= 
onderstellen, dat ín dat geval zich als gevolg van (bev.) de 
tearing=mode magnetische eilanden ontwikkelen die dezelfde dike 
te hebben als de schijf, en zo magnetische velden genereren, 
die hetzelfde zijn als turbulente-velden met de golflengte 4he 
Omdat de maanetische turbulentie ín de radiële richting niet 
direct wordt beïnvloed, veronderstellen we, dat ook de viscosie= 
teit (zie vol.(7)) niet wordt veranderd, 


Zo vinden wet: 


de b 
s 4 a 
De (B, we 5) (15) 


Deze vergelijking lossen we als volgt op: 


1) We zoeken een oplossing voorde homogene vergelijking 


du 4 
ge = 0073 7 U, (16) 
dat is U = C Ede (17) 


waarin C een integratie constante is, 
2) De oplossing van de inhomogene verbelijking wordt dan gege= 


ven door: 


B, = Cr0°75 4 0,73 { E 2 dr. (18) 
s de Se 4 
r : bad 


o 


Uit de randvoorwaarde, B (r) = 0, volgt direct dat C nul moet 


zijn, Zo vinden we. uiteindelijk voor het verticale magnetische 


veld: 
vn et te at 
z= 75 s 3 \T 573) " H073T Ee 
r B Ü oa 
(19) 
=S (1,196 Ne 
T 





4á 


waarin y = r/r, de dimensisloze straal:is, 


10 








1 4.2 4.4 1 5 wr 
Figuur 17 De radiële afhankelijkheid van B. 


De grens van het overgangsgebied is gedefiniëerd als het eerste 
nul=punt van vergelijking (19) en een numerieke benadering voor 


dit nulpunt is (zie ook figuur 17): 


r, = 1,624 re (20) 


De vergelijkingen IIll=(7), III=(13) en III=(14) kunnen nu 
opnieuw numeriek worden opgelost, Omdat het overgangsgebied 
klein is, moeten we in vergelijking IIIl=(7) nu wel de gradiënt 
van de magnetische druk meenemen, Er is dan sprake van een stele 
sel partiële differentiaal vergelijkingen, In de volgende para= 
graaf beschrijven we hoe we dit stelsel opgelost hebben en zul= 


len we de resultaten bespreken, 
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IV=2 De numerieke berekening en de resultaten, 


We moeten nu het volgende stelsel vergelijkingen oplossen 
(c.f . pag.30 t/m 33): 
2 2 


v B 2 
ae 2 EN dte op DES 
mes 25e + est = 5 + A dd DIE) IIIl=(7) 


le 


z B 
e(z) = e(z=0) = er | (z ae + zie) dz en _III=(13) 
o 


áT 
r 
z ATO ep Ì 
2) = \ En - vib) ez. III-(14) 


De temperatuur en de dichtheid in het equatorvlak worden gege= 
ven door de oplossing van de Ichimaru=schijf ín vergelijking 
(11); Het effectieve geleidingsvermogen is gegeven in vergelij- 
king (6), L = 4h (cf page43); B, volgt uit vergelijking (19); 
En \_ volgt uit de Ichimaru-schijf oplossing en de wet van be= 
houd van massa (zie pag,10 de tweede veronderstelling). 

Zoals we in de vorige paragraaf al opmerkten zorgt de ale 
lerlaatste term in vergelijking IIlI=(7) voor problemen: deze 
maakt dat het een stelsel partiële dífferentiaalvergelijkingen 
is, Dezelfde methode van oplossen als in paragraaf III-1 is nu 
niet direct mogelijk, omdat de radiële afgeleide van de magne- 
tischedruk nog onbekend is, 

We lossen dit probleem op. door het stelsel vergelijkingen 
simultaan voor drie dichtbijelkaar gelegen waarden van r op te 
lossen; dat zijn r =r= ér, Tr en EE + re We gebruiken nu 


de Taylor=reeksontwikkeling voor functies van r‚ dei, 
+ + 2 
f(r ee Sr) = f(r) em ârf'(r) + +$6r Pt'(r) ee (21) 


evenzo f'(r pd Sr) = f'(r) E Erft!(r) …… (22) 
Als we nu f(r — &r), f(r) en f(r + 8r) kennen, kunnen we met 

de vergelijkingen (21) en (22) ook f'(r — êr), f'(r) en 

f'(r +8r) bepalen. We vinden dan tot op eerste orde in r: 


f'(r Ss &r) = áf(r) _ sf(r Kea êr) _ f(r + ôr) , (23a) 


2ór 


f'(r) ne fr + AT In r=êr (23b) 


f'(r + 6r) == 3f(r + ôr) + PLE mcl 8r) ams 4f(r) (23c) 


2ör 
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Voor &r kiezen we waarden van de orde van 107Sr; bij het ge= 
bruiken van veel kleinere waarden voor Sr moeten we voorzichtig 
zijn, omdat de verschillen in vergelijking (23) dan van de orde 


van of kleiner dan de machine=nauwkeurigheid kunnen worden (zie 





Stoer, 1972, pag.á e,Ve). Op deze manier kunnen we, als we het- 
magnetische veld op deze drie plaatsen kennen, de radiële afge- 


leiden bepalen en de radiële drukgradiënt berekenen, De verge | 


lijkingen III=(7), III=(13) en IIl=(14) kunnen we dan numeriek © 
integreren, 
De numerieke berekeningen hebben we weer met het IBM pro= ° 





gramma ESMP III op de DEC-20 van het CIW gedaan (zie Appendix 
A). Het specifieke CSMP-programma voor de boven beschreven me 
thode geven we in Appendix C, De dimensieloze eenheden, die we 
in dit programma gebruiken zijn op de gebruikelijke manier van 
de voorgaande paragrafen gedefiniëerd, behalve in de verticale 
coördinaat, die gedefiniëerd wordt als z/h, waarin h nu niet 
de halve dikte van de schijf is, maar de schaalhoogte van de 
druk (eef. pag.33/34): 


Jor Ì 
= 72E RT 
hz ECT Ie Ve III=(16) 


Net als in de berekeningen uit paragraaf III=1 moesten we 
ook nu de parameter C gebruiken om het radiële massatransport 
ongeveer gelijk aan n te houden, zie ook pagina 34, 

Omdat we speciaal geïnteresseerd zijn in de toepassing van 
het model op Her X=1, hebben we de berekeningen gedaan met 


waarden die typisch zijn voor dít systeem: 


M= 1,4 M (m = 1); n= 1007 a/s (fh = 1); 


3) 
en P= 1,24 3, 


Op pagina 28 lieten we al zien, dat hieruit ook een relatie 
tussen het magnetische dipoolmoment en de fastness kan worden 
afgeleid, dei. @ 


« 


k= 1,7 2030 0/6 em. Ill=(2) 
8 
De resultaten, 


In figuur 18 geven we twee typische voorbeelden van de ge= 
vonden verticale structuren; beide zijn berekend bij de fast= 
ness 0,9, wat voor het dipoolmoment een waarde van 1,5 10% em® 
geeft, Figuur 18a) is berekend voor r = 1,06 r_‚ d‚i, binnen de 
corotatiestraal, en 18b) voor r:= 1,08 Fos dede buiten de coro= 


tatiestraal, 
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Figuur 18 De verticale structuur, a) r<r 
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De parameter C kon in onze berekeningen steeds kleiner dan 
10 gekozen worden, zodat de schijf dus niet zo extreem gecome 
primeerd wordt als ín het model van Ghosh en Lamb (zie pag.35). 


Een reeks berekeningen als in figuur 18 voor waarden van r 





die het hele overgangsgebied (rs r<r= 1,624 r) bestrijken, 
geeft ons in het overgangsgebied ke In figuur 19 schetsen we de 
zo verkregen resultaten voor vier verschillende waarde van de 
fastness, mende. ah we = 0.6, 0e7, 0,8 en 0,9, Volgens pe 
vergelijking III=(2) correspondeert dit voor Her Xe1 met magne= 
tische dipoolmomenten van respectievelijk 0,94, 1,1, 1,3 en 


1,5 10°? G em°. 








ee 
co 
1.2 : 44 — r/r, en, 8 






4 E 
1.2 ad 1.4 rs r/m 1é 


Dad 


b) 
Figuur 19 a) en b) De resultaten voor ze in het overgangsgebied, 
Zie ook de volgende pagina, 





á9 





r 
rn 1.é $ 


10 





d) B Yo an 14 > tg Mb 5 


Figuur 19 c) en d) De resultaten voor y& ín het overgangsgebied, 
Zie ook de vorige pagina en merk op, dat binnen de co= 
rotatiestraal Yp maximaal ongeveer 50 wordt en dan 
verder naar binnen constant blijft, 

- Vervolgens kunnen we met vergelijkingen Il-(45) en II=(46) 


het totale accretiekoppel op de neutronenster berekenen: 
° T . 
K= M(Gmr_)? + { ° le), (rr? or (24) 
o É $ z 
ao 


dar r 


De resultaten van deze berekening Geven we ‘fh fiduür 20, We 
schatten de onzekerheid in deze resultaten op ongeveer 10%, dit 


voor zover ze veroorzaakt ís door de numerieke berekening. 
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ob og ï 


Figuur 20 De berekende accretiekoppels, 
De lijn is een op het oog getrokken beste lineaire fit, 
Zo vinden we voor Her X=1, waarop het koppel volgens de 
waarnemingen nul of bijna nul moet zijn, een fastness 
van ongeveer 0,87, 
In figuur 20 zien we, dat de fastness van Her X=1 ín dit model 
ongeveer 0,87 moet zijn, We vinden zo een redelijke overeen 


stemming tussen het model en de waarnemingen voor: 


Kd 
M= 1.4 Ms M= 1017 o/s; 


0’ 
5, = 0,87 en Me 1,5 10% G em”, 


Een belangrijk bezwaar tegen deze berekeningen is helaas, 
dat de magnetische spanning ten gevolge van de interactie van 
het magnetische veld van de neutronenster met de accretieschijf 
erg groot wordt Ge is vrij groot): de invloed hiervan op de 
radiële structuur van de schijf is dan niet meer te verwaarlo= 
Zen, zodat onze resultaten niet consistent zijn met de gedane 
veronderstellingen, 

Wel laten deze berekeningen opnieuw zien, dat er een krie 
tische fastness kan bestaan waarvoor het totale accretiekoppel 


op de neutronenster hul wordt, 
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Va1 Het model van Anzer en Börner. 


Nu we het model van Ghosh en Lamb uitvoerig hebben bestu= 
desrd, willen we nog twee andere modellen bespreken, We begin= 
nen met het model, dat door Anzer en Börner, 1980 en 1983, is 
voorgesteld, Zij veronderstellen, dat het magnetische dipoolmoe 
ment van de neutronenster niet parallel ís aan, maar loodrecht 
op de rotatieas van de neutronenster en de accretieschijf staat, 
Verder veronderstellen zij, dat de accretieschijf perfect diaa 
magnetisch is, Dat houdt in dat het magnetische veld van de 
neutronenster niet in de accretieschijf kan doordringen en dus 
ook. dat het geleidingsvermogen van de accretieschijf oneindig 


is. In figuur 21 is deze configuratie geschetst, 





… Fig. 21 Model of the magnetosphere-disk configuration. The 
upper part shows the side view of the magnetic dipole field, the 
lower part gives the projection of the field lines when the dipole à 
is seen end-on. (Anzer À Börner, wao 


Wel bestaat er natuurlijk een snelheidsverschil tussen de 
Keplers roterende accretieschijf en de starre magnetosfeer; dit 
snelheidsverschil op de grens van schijf en magnetosfeer ver= 
oorzaakt een Kelvin=Helmholtz instabiliteit, Door de Kelvin= 
Helmholtz instabiliteit worden de materie uit de schijf en het | 


magnetische veld vermengd, Als dit binnen de corotatiestraal 
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gebeurt, zal de materie via de magnetische veldlijnen op de 
magnetische polen van de neutronenster vallen, Als dit buiten 
de corotatiestraal gebeurt, zal de materie door het magnetische 
veld worden weggeslingerd, 

In het eerste artikel (Anzer & Börner, 1980) bestuderen 
Anzer en Börger de Kelvin-Helmholtz instabiliteit tussen twee 
media met dezelfde geluidssnelheid, Ce, maar verschillende 
dichtheden, e4 in de schijf en e‚ in de magnetosfeer, Zij vin= 
den, dat voor grotere snelheidsverschillen tussen beide medía 
de Kelvin=Helmholtz instabiliteit wordt onderdrukt; dit gebeurt 
als 

1 
Av Cosx >} 2E) e« (1) 


Hierin is « de hoek tussen het snelheidsverschil en de golf vec= 


tor van de verstoringe 


In het tweede artikel (Anzer & Börner, 1983) breiden Anzer 


en Börner dít uit door ook rekening te houden met een magne- 
tisch veld ín het tweede medium (tangentieel aan het grensop= 
pervlak tussen beide media), Uit de eis voor druk evenwicht 
volgt dan, dat de magnetische druk boven de schijf gelijk moet 
zijn aan de gasdruk in de bovenrand van de accretieschijf, Ook 
dit magnetische veld kan de Kelvin=Helmhotz instabiliteit stan 
biliseren, Voorwaarde voor stabiliteit is nu (zie Anzer & Bör- 
ner, 1983): 


av cosayn/2 Vage (2) 
Hierin ís Va de Alfén=snelheid ín de magnetosfeer: 
B 
Va 7 ne,” (3) 


De eis van drukevenwicht kunnen we schrijven als 
2 2 
Ots = ZOVt n (4) 


Met vergelijking (4) kunnen we de voorwaarde voor stabiliteit 
uit vergelijking (2) herschrijven en we vinden dan dezelfde 
voorwaarde voor stabiliteit als ín het geval zonder magnetisch 
veld, dei, vergelijking (1)e 

Bovendien zullen de magnetische veldlijnen sterk verstoord 
worden als Av> c,, omdat dan i 


2 
ea)? > 5. (5) 
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Dit leidt tot een sterke opwinding en dus tot grote azímuthale 
componenten van het magnetische veld, díe de Kelvin=Helmholtz 
instabiliteit nog eens extra onderdrukken, 

Bovendien volgt al uit de voorwaarde in vergelijking (1), 
dat voor een typisch dichtheidsverschil, p‚ = 100 @,‚ alleen 
een relatief smalle ring rond de corotatiestraal instabiel is, 
Hierdoor is het accret{ekappel, dat het gevolg is van de boven 
beschreven interactie tussen de schijf en de magnetosfeer, in 
eerste benadering eehvoudig te berekenens 
=de inval van de materie die binnen de corotatiestraal met het 
magnetische veld wordt gemengd, geeft een positieve bijdrage 
aan het totale koppel: 

e 2 : 

Kace * Pinte” (6) 
Hierin is Je de totale invallende massa per seconde en is ® 
de rotatie snelheid van de neutronenster en dus ook van de ma= 
terie, omdat deze immers vlak bij de corotatiestraal uit de 
schijf wordt gehaald, 
„de materie die buiten de corotatiestraal uit de schijf wordt 
gehaald.en wordt gedwongen tot corotatie met de magnetosfeer, 
voelt een effectieve zwaartekrachtsversnelling, die naar buiten 
is gericht, Hierdoor wordt deze materie dus naar buiten ge= 
transporteerd en versneld; op een gegeven moment heeft deze ma= 
terie zoveel energie gewonnen, dat ze niet langer aan de magne= 
tische veldlijnen gebonden kan worden, Als dit gebeurt op de 


straal Feg Een = ) wordt de impulswinst van deze materie 


IE 
eq co 
en dus het koppel, dat de rotatie van de neutronenster breekt, 


gegeven doors 


K =n Rr (tex) (7) 
b out “co eq” 
. 
Hierin ís Mit het totale massaverlies van de schijf buiten de 
corotatiestraal, à 
Ld Kd 
Anzer en Börner veronderstellen nu, dat M en M aan 
e in aut 


elkaar gelijk zijn en noemen dit M, Verder veronderstellen zij, 
dat er elders in de accretieschijf geen materie meer de schijf 
kan verlaten, ook niet aan de binnenkant van de accretieschijf, 
Sommeren van vergelijkingen (6) en (7) geeft dan het totale 


accretiekoppel op de neutronensters 


Kk = Mar 2 (2 - Xen) _ (8) 





BÁ 


Anzer en Börner zijn niet in staat om Xag te berekenen en 


definiëren daarom de vrije parameter sz 


Ze 2 
xs 2e Kea* | (9) 
Uit vergelijking (8) voor het totale koppel op de neutronenster 
kunnen we nu de periodeafname van de neutronenster bepalen (mats 
pageá): 


Pp =.8,9 10° ars, m2/3 â p7/3 s/ jr. (10) 


Anzer en Börner veronderstellen nu, dat de Röntgenlichtkracht 
gegeven wordt door 0,1 me? voor de meeste neutronensterren 

komt dit overeen met de vrijgekomen potentiële energie van de 
invallende materie. Vergelijking (10) kunnen we dan schrijven 


als: 


=P = 8,9 10° « ien n2/3 L_ P7/3 o/sr. (11) 


37 


Kd 
In figuur 22 zijn de P's van een aantal Röntgenpulsars uitgezet 


V4LREN De rechte lijn ín deze figuur geeft de theo= 


retische relatie volgens vergelijking (11) meto{= 0,05 en Ius = 


tegen de P 


m = 1, De meeste waarnemingen komen dus redelijk overeen met de 


theoretische voorspelde relatie (zie Anzer & Börner, 1983). 


tog (-P) 






Cen X-3 
SMC X-1 + 


tog! p/3 L37) « 


2 n hk 
Fig à A plot of P versus aP”/*M (with a=0.05) compared to the 
observations. Vela X-1 is plotted in its spin-up phase , 
The dotted line corresponds to a white dwarf of 1M 


Er zijn echter een aantal belangríjke problemen in dit 


model: 
1) Zoals Anzer en Börner, 1983, zelf al opmerken, is de overeen 


komst tussen de voorspellingen van hun model en de waarne= 


mingen niet zo verwonderlijk, omdat Xeq (en «) neg een vrije 
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parameter is, zodat ín feite met de juiste keuze voor “eg bijna 


alles verklaard kan worden. 


2) In dít model ís aangenomen, dat al het massaverlies uit de 
schijf in de ring rond de corotaiestraal plaatsvindt, Dat 
betekent, dat er in de schijf binnen de corotatiestraal geen 
massatransport meer mag plaatsvinden óf dat binnen de ring geen 
schijf meer beëtaat, meaewe de binnenrand van de schijf wordt 
nier gegeven door de magnetosferische straal maar door de bin= 
nenrand van de ring waarin het massaverlies vit:de schijf 
plaatsvindt, In elk geval zal dit grote massaverlies de radiële 
structuur van de accretieschijf sterk veranderen, bijvoorbeeld 


zoals in de grenslaag van het model van Ghosh en Lamb, 


3) Van groot belang ín dit model is, dat de accretieschijf in= 
staat ís zich volledig van het magnetische veld van de neu= 
tronenster af te schermen, Dit betekent onder andere, dat het 
geleidingsvermogen oneindig (of heel groot) moet zijn, Als het 
geleidingsvermogen gegeven wordt door zijn klassieke waarde, Ís 
dit een redelijke veronderstelling, Maar we verwachten niet, 
dat het effectieve geleidingsvermogen van de schijf inderdaad 
door zijn klassieke waarde gegeven wordt, In een accretieschijf 
moet de viscositeit immers vele orden groter zijn dan zijn 
klassieke waarde en verwachten we, dat de processen die deze a= 
nomale viscositeit veroorzaken, ook een anomaal geleidingsver= 
mogen geven (c‚f‚ page 35 en 38). In het schijf model van 
Ichimaru, 1977 (zie paragraaf IV-1), wordt zo inderdaad een ef- 
fectief geleidingsvermogen gevonden, dat vele orden kleiner is 


dan het klassieke geleidingsvermogen. 


4) De accretieschijf kan zich níet binnen de straal waar de 
magnetische druk in de magnetosfeer gelijk wordt aan ce gas= 
druk in het centrum van de schijf, uitstrekken, In het model 
van Anzer en Börner (in feite in elk model, waarin accretie op 
de neutronenster plaatsvindt) moet de corotatiestraal dus bui 
ten deze straal liggen, Geheel analoog aan de afleiding van de 
relatie tussen de fastness en het magnetische dipoolmoment van 
Her Xe1 (zie pag,28) kunnen we uit deze voorwaarde een. bovene 


grens voor het magnetische dipoolmoment van Her X-1 vinden; 


HK < 1,9 10°° Gem®, : (12) 


Dit is duidelijk kleiner dan het magnetische dipoolmoment, zoe 
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als dít volgt uit de waarneming van de ciíiclotron=emissielijn in 
het harde Röntgenspectrum van Her X=1.(zie pag.27/28), 

N.B. dít probleem is niet specifiek voor dit model, maar geldt 

voor elk model, dat de accretie op Her X=1 met een schijf ver= 


klaart! 
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V-2 De scheve rotator ('oblique rotator'),. 


In alle voorgaande paragrafen hebben we steeds het prom 
bleem in een speciale geometrie bestudeerd, waarin het magnetie 
sche dipoolmoment en de rotatieassen van de neutronenster en de 
accretieschijf parallel waren of het magnetische dipoolmoment 
loodrecht op de beide rotatieassen, die weer parallel waren, 
stond, We willen nu de meer algemene situatie bestuderen, waar= 
in het magnetische dipoolmoment en de beide rotatieassen wille= 
keurige hoeken met elkaar maken, Daartoe beschouwen we eerst 
het probleem van een niet roterende neutronenster, waarvan het 
magnetische dipoolmoment een hoek X maakt met de rotatieas van 
de accretieschijf (zie figuur 23), die aan de binnenkant be= 


grensd ís door de magnetosferische straal, En Lipunov, 1978 





Figuur 23 De geometrische configuratie, De schijf ligt in het 
XY=vlak en ú is het magnetische dipoolmoment van de 


neutronenster, 


(zie ook Kundt & Robnik, 1980), berekende de magnetische velden 
voor het twee=dimensionale analogon van dit probleem met een 
oneindig dunne schijf met oneindig geleidingsvermogen, De vol-= 
ledige drie=dimensionale oplossing wordt gegeven door Aly, 
1980, opnieuw voor een oneindig dunne en oneindig goed gelei 
dende schijf en déze oplossing verschilt nauwelijks van die van 
Lipunov, 1978, voor het tweedimensionale analogon. 

We zullen hier níet ingaan op de berekeningen en geven 
alleen de belangrijkste resultaten, Het totale magnetische veld 
(dipoolveld + veld van de in de schijf geïnduceerde stromen) is 
volgens Aly, 1980 (in sferische coördinaten r, 9 en @): 
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B, „en [stnocos(p-ptsinX + 2eoskleoslatanx + 0) ‚ (13a) 
r 
Bg = Ee [-eoeseostg-a ear + 2 2 (13b) 
r 4 2 2 rs+r 
ri 2eoskeind „tanixi mn Seren B, YX. a — 1)) 
cos 8 r 
m 
en 
m ‘ 
B == sinlp=p')sink, (13e) 
P 5 
Hierin iss 
2 } 
X= Er, cos8 ECE CPR en (13d) 
r-nt Er Zy 
Y= rr? zele + 4r2rcos”0) . (13e) 


Figuur 2á geeft een globaal beeld ven dit magnetische veld, 


VD ee 
7 is Wid EN El d 





Figuur 24 De globale structuur van het magnetische veld volgens 


_ Dams 
Aly, 1980, b) in het vlak door M en Sg, 


De karakteristieken van deze oplossing zijn (zie ook Kundt & 
Robnik, 1980): 
1)-er zijn twee neutrale punten, N° en N , op het oppervlak van 
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de accretieschijf en de coördinaten van deze punten zijn: 


0%; 
… d 
N r idem. p= en z 0 


= 
H 
H 


r (1 + al cotan2x)?, p= 0 enz 
k n (14) 


2) voor F&r, vinden we het ongestoorde dipoolveld terug en 
voor rr, wordt het veld dat met 4, ís geassociëerd (de 
termen met cos in vergelijking (13)) gereduceerd tot een qua=- 
drupoolveld, omdat de geïnduceerde stromen in de schijf een díi-= 


poolmoment hebben, dat precies tegenovergesteld is aan Hye 


3) de magnetische druk onder en boven de accretieschijf is niet 
overal hetzelfde en dit drukverschil resulteert in een kope 


pel op de accretieschijf, dat gegeven wordt door (zie Aly,‚1980) 


KA 
3n 


EN 
K = 


ces X Aoxn)s (15) 


JP 
3 N[ N 


waarin f, en Au de eenheidsvectoren zijn in de richting van 
respectievelijk het magnetische dipoolmoment en de rotatieas 


van de accretieschijf, 


Hoewel deze berekening gedaan ís voor een niet-roterende 
neutronenster, veronderstellen we nu, dat de oplossing níet 
zal veranderen als de neutronenster wel roteert, omdat het ge= 
gied, dat we bestuderen meestal ver binnen de lichtcylinder zal 
liggen. Het koppel uit vergelijking (15) geeft ons dan op elk 
moment het op de schijf vitgeöefende koppel, dat nu dus perio= 
diek variöert, De meest algemene geometrie is nu in figuur 25 


geschetst, De rotatieas van de accretieschijf maakt een hoek 





Figuur 25 De algemene geometrie, 
Ag Ne en Au zijn de eenheidsvectoren in de richting van respe 
de rotatieas van de schijf en de neutronenster en het dipoolmo= 


ment 
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met de rotatieas van de neutronenster en het magnetische dipool 
moment maakt een hoek met de rotatieas van de neutronenster, 
Het over één rotatieperiode gemiddeldes koppel op de schijf ís 
nu (zie Lipunov & Shakura, 1980): 








2 
WK = 5 22 cosa(Icos fs = 1) (Axhg)e (16) 
m 
Lipunov en Shakura, 1980, geven ook de gemiddelde energie van ® 
de magnetische dipool van de neutronenster in het magnetische 
veld van de in de schijf geïnduceerde stromen: é | 
âme 2 2 2 
KE) = ze Eel 2eos fp + (1 - Jcos“B)sin ox). (17) 
37 "2 
m 


De stabiele positie van de schijf wordt nu gegeven door de po= 
sitie waar <E> minimaal is. Volgens vergelijking (17) is dat: 


«= 0 als (3) arccos atd of 
: (18) 
ï - ‘ 
«== 5 als (3 <arccos 3 *, 


Het totale koppel op de neutronenster is nu de som van het ma= 
teriële koppel als gevolg van de koppeling van de materie aan 
het veld op de binnenrand van de schijf en het magnetische kop= 
pel, dat volgt uit het op de schijf uit geöefende koppel (actie 
= ereactie). Over het algemeen zal het magnetische koppel veel 
groter zijn dan het materiële koppel en als we het koppel dan 
splitsen in een bijdrage evenwijdig met de rotatieas van de 


neutronenster en een bijdrage loedrecht op deze as, vinden we: 


Kd 
Ki, =/GMr cosa M en 


> á 
7 


(19) 
cosol Jcos - 1) (Ax). 


nF 
BNI N B 


K, geeft nu de spinup (of spindown) van de neutronenster en K 
doet de neutronenstsr precederen, Zo vinden we de volgende re= 
sultaten (zie o.a, Lipunov & Shakura, 1980): 

1) als P) arccos 3 °, zal S&O en is Ki, = /Mr, M, maar als 
(3< arccos 3=Â, heeft de schijf zijn stabiele positie voor 
x= Ee en vinden we een veel kleinere K,, of zelfs een negatieve 
Kyje 
Her Xe1 en Cen X=3, Nauwkeuriger berekeningen van Lipunov, 


Dit verklaart de ongewoon langzame perigdeaf name van o,a, 


Semnev en Shakura, 1981, laten echter zien, dat onafhankelijk 
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van f het binnenste deel van de accret{ieschijf naar = 0 nijgt 


als { verder naar buiten niet te groot wordt, 


2) de precessie periode van de neutronenster wordt gegeven door 
(zie ook Líipunov & Shakura, 1982): 


3 1 2 _1 2 3 o=eÎ 
Tr = 6,2 10” sin 2% (3cos B — 1) F3o Tas PaP jre (20) 


Hierin is Ra = r,/1ofem, We zien nu, dat deze precessie periode 
voor een langzaam roterende neutronenster met een sterk magne= 
tisch veld van de orde van enkele weken tot enkele dagen kan 


zijn (bev, Vela X=1), 


Hoewel in dit model opnieuw is verondersteld, dat de 
schijf een oneidig geleidingsvermogen heeft, verwachten we niet 
dat de globale structuur van het magneetveld sterk zal verande= 
ren, als deze veronderstelling niet geheel juist zou zijn, om= 
dat in het model van Ghosh en Lamb ook maar een klein deel van 
de magnetische flux in de schijf kan doordringen, We verwachten 
dat ook in een dregelijk model de boven vermelde resultaten 
geldig blijven, 

Bovendien is het mogelijk in het kader van dit model vele 
eigenaardigheden van de Röntgenpulsars te verklaren (zie o.a, 
Lipunov & Shakura, 1980) en vindt ook een object als SS 433 
een natuurlijke verklaring binnen dit model (zie Lipunov & 
Shakura, 1982), Aly, 1980, geeft ook een qualitatieve beschrije 
ving van de accretiestroming uit de schijf naar de neutronen= 


ster, hoewel dit probleem verre van definitief opgelost is, 
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VI Slotopmerkingen en aanbevelingen voor toekomstig onderzoek, 


In de voorgaande hoofdstukken hebben we laten zien, dat de 
modellen van Ghosh en Lamb en van Anzer en Börner nog een groot 


aantal problemen oproepen: 





Het grootste probleem ín het model van Ghosh en Lamb ís de 
schending van het verticale hydrostatische evenwicht, Deze . 
hangt ook nauw samen met de omstreden afleiding van het effec= 
tieve geleidikgsvermogen, De door ons voorgestelde veranderin= 


gen hebben deze twee bezwaren wel kunnen opheffen, maar zijn op 





een ander punt zelf inconsistent gebleken, Verder onderzoek en 
met name het tegelijk oplossen van de radiële, azimuthale bewe= 
gingsvergelijking en de vergelijking voor het hydrostatische 
evenwicht is nog nodig, Ook betere afleidingen van het effec 
tieve geleidingsvermogen zijn nodig, waarbij speciale aandacht 
aan de relatie tussen dit effectieve geleidingsvermogen en de 
anomale viscositeit van de accretieschijf geschonken moet wor= 
den. Zolang deze relatie nog onbekend is, kan het zinvol zijn 
te bedenken, dat de electrische resistieviteit en de kinemati= 
sche viscositeit dezelfde dimensie hebben, zodat een eenvoudige 


evenredigheid tussen beide voor de hand ligt, 


In het model van Anzer en Börner ís een goede berekening 
van Ken meer dan nodig en meot ook de radiële structuur van de 
betrokken ring rond de corotatiegstraal worden bestudeerd, Zo= 
lang dit nog niet is gebeurt, is het in feite niet mogelijk 
veel zinnigs te zeggen over de mogelijkheden of onmogelijkheden 


van dit model, 


Bovendien lijkt het onwaarschijnlijk, dat één van deze 
twee modellen alle waargenomen periodeaf names van de Röntgen= 
pulsars kan verklaren, Beide gaan immers uit van een heel spe= 
ciale geometrie, waarin het magnetische dipoolmoment parallel 
is aan of loodrecht staat op de rotatieas van de „neutronenster, d 
Waarnemingen van de hoek tussen dipoolmoment en rotatieas bij 
‘gewone sterren! geven een verdeling met twee brede pieken (on= 
geveer 20°) rond de 0° en de 90°, die ongeveer even hoog zijn 
(zie Preston, 1971). In eerste instantie verwacht je dus ook 
een dergelijke verdeling voor de Röntgenpulsars, Bovendien zijn 
in de beide modellen de rotatieas van de schijf en die van de 


neutronenster parallel en ook dit lijkt niet algemeen noodzake= 
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lijk te zijn voor Röntgenpulsars. 

Het model van de scheve rotator, waarin alle drie de assen 
vrij gekozen kunnen worden, verdient dan ook de voorkeur, Bo 
vendien is het zeer waarschijnlijk, dat de modellen van Ghosh 
en Lamb en van Anzer en Börner, als ze uiteindelijk volledig 
consistent zijn, eenvoudig als limiet gevallen in dit veel ale 
gemenere model kunnen worden ingepast, 

Een belangrijk probleem blijft wel het massaverlies uit de 
schijf en de daarop volgende stroming naar de neutronenster, 
Verder onderzoek langs de door Aly, 1980, aangegeven weg is dan 
ook hard nodig, 

Ook mogen we hier het twee=stromen-accretiemodel van Lipu= 
nov, 1982 a en b, niet onvermeld laten, In dit model bestudeert 
Lipunov de gecombineerde accretie door Roche-lob-overstroming 
(schijfaccretie) en door een sterrewind (sferische accretie), 
Door de combinatie van deze twee accrtiestromingen kunnen vele 
facetten van het onregelmatige periodegedrag van de Röntgenpul= 
sars worden verklaard, 

Tot slot nog een opmerking over de waargenomen electronen= 
eyclotron=emissielijn in het harde Röntgenspectrum van Her Xe: 
uit een alternatieve verklaring voor het ontstaan van deze 

lijn van Bonazzola, Heyvaerts en Puget, 1979, kan een dipool 


za e worden afgeleidj een resültaat, dat in 


moment van 1,á 10 Gem 
tegenstelling tot het resultaat van Trümper et al, 1978, goed 
is in te passen in alle drie de behandelde modellen (cef. pagi= 
na 28/29 en 55/56), Een definitieve bevestiging van één van 


beide resultaten of een ander resultaat is zeer gewenst, 
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Appendix Ae 


CSMP III is de drede versie van een IBM=-programma voor het 
simuleren van continue systemen, In deze bespreking van CSMP 
zullen we ons beperken tot die details, die nodig zijn om de in 
de volgende paragrafen te bespreken CSMP-programma's te begrij= 
pen. Meer details kunnen gevonden worden in de IBM publicatie 
no. SH19=7001=3 “Continuous System Modelling Program Program 
Reference Manual", 

Met CSMP is het mogelijk een stelsel niet-lineaire, gewone 
differentiaal-vergelijkingen op te lossen, Hiervoor zijn een 
achttal verschillende integratisroutines beschikbaar, Voorzover 
de rekentijd dit toeliet, hebben we steeds gebruikgemaakt van 
de vierde orde Runge-Kutta methode met een variabele integra= 
tiestap; in de andere gevallen, als deze methode teveel tijd 
kostte, hebben we dezelfde mathode met een vaste integratiestap 
of de gewone Euler methode gebruikt, De details van deze metho= 
des en hun definities zijn te vinden in bijvoorbeeld Stoer, 
1973 Hoofdstuk 7. 

„Een CSMP programma bestaat uit een aantal CSMP statements 
en eventueel een aantal gewone FORTRAN statements, Deze state 
ments worden door CSMP gesorteerd en daarna vertaalt in een 
FORTRAN programma, UPDATE, waarna dit programma door de FORTRAN 
compiler in machinetaal kan worden vertaald en vervolgens kan 
worden uitgevoerd, i 

We kunnen een CSMP programma opdelen in drie delen: 
INITIALe= Een eerste deel waarin de begincondities en een aantal 

parameters kunnen worden vastgelegd, 

DYNAMIC Dit deel vormt de kern van het programma, die aangeeft 
hoe de diverse grootheden samenhangen en welke groot= 
heden geïntegreerd moeten worden, 

TERMINAL-In dit deel worden de eindcondities vastgelegd (duz, 
wanneer en onder welke voorwaarden de integratie ge= 
stopt moet worden) en kan worden aangegeven van welke 
grootheden men de resultaten wil zien, De resultaten 
kunnen in de vorm van een tabel of een plot worden af= 
gedrukt, 

De opdeling van het CSMP programma ín deze drie delen kan ( 

hoeft niet) aangegeven worden door de statements INITIAL, 

DYNAMIC of TERMINAL aan het begin van de respectievelijke delen 


te plaatsen, 





A2 


Naast de gewone FORTRAN statements zullen we de volgende 


typische CSMP statements gebruiken: 


PARAMETER= Dit statement dient am aan te geven, dat er na een 
aantal parametrs worden gedefiniëerd, Bijvoorbeeld 
A en B zijn gelijk aan 1: 
PARAMETER A=1., B=1. 





INTGRL Dit statement dient om aan te geven, dat een groote 





heid (bev, A) berekend moet worden door integratie v 
van een andere grootheid (bev. B) bij een vaste bee 
ginvoorwaarde (bev, A(t=0) = Af): ii 
A=INTGRL(Af,B) 
NOT Dit statement kan men gebruiken om te constateren, 
dat een grootheid van teken verandert: 
Y=NOT (XxX). 
betekent, dat Y 0 is, als X groter ís dan 0, en Y 1 
is, als X kleiner is dan 0, 
TIMER= ‘Met het TIMER statement kunnen de integratjestap= 
grootte, DELT, de maximale 'tijd' waarover geïnte | 
greerd wordt, FINTIM, en de stapgrootte waarin de 
uitvoer moet geschieden, OUTDEL of PRDEL, worden 
vastgelegd, Bijvoorbeeld: 
TIMER DELT=0,01, FINTIM=1., OUTDEL=O 1 
METHOD Voor de definitie van de te gebruiken integratie= 
methode: 
METHOD RKS = gebruik Runge=Kutta met variabele stapgrootte, 
METHOD RKSF= idem met vastte stapgrootte, 
METHOD RECT= gebruik Euler methode, 
FINISHe Met dít statement kan worden aangegeven, dat de ine 
tegratie onder bepaalde voorwaarden gestopt moet 
worden, Bijvoorbeeld als de integratie gestopt moet 
worden als B van teken veranderts ú 5 
FINISH B=0, 8 
OUTPUT. Hierna kan men de grootheden aangeven, waarvan men 


een plot of een tabel wil hebben, 


INTERACTIVE-Dít statement, dat direct na het laatste PARAMETER 
statement geplaatst dient te worden, geeft weer, dat 
men het programme in de interactieve mode wil ge= 


bruiken, 
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In de interactieve mode kan men het CSMP programma vele 
malen na elkaar uitvoeren, zonder dat het CSMP programma steeds 
opnieuw vertaalt moet worden, Voor elke 'run' kan men dan de 
diverse parameters definieren, zodat de interactieve mode erg- 
geschikt is om een zelfde set vergelijkingen met steeds veran-. 
derende parameters of begincondities snel door te rekenen, 

Van elke set differentiaalvergelijkingen kunnen we eenvoum 
dig een blokschema maken, waarbij elk blok in een CSMP of 
FORTRAN statement kan worden vertaalt, Dit ziet er voor de set, 


Ede cosa/B _gmet randvoorwaarde A(lt=0) = 0, 
dB A 
gE smet randvoorwaarde B(t=0) = 0, 


dan als volgt uits 
Het blokschena: 





Het CSMP programma: 


INITIAL 
PARAMETER Ad=f, Bf=d 

DYNAMIC 

DA=COS(SQRT(B)) 

DB=EXP(A) 

A=INTGRL(A#,DA) 

B=INTGRL(8f,DB) 

TERMINAL 

TIMER DELT=d.f1, FINTIM=1e, OUTDEL=d. 1 
OUTPUT A,„B 

END ° 
STOP 
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Appendix Be 


In deze appendix bespreken we het het CSMP programma, 
waarmee we de verticale structuur van de accretieschijf berske= 
nen op de manier zoals die op pagina 30 t/m 33 is gegeven, Dat 


betekent, dat we de volgende vergelijkingen moeten oplossen: 











2 
B 
2 2 -? GM & 
vs + 25T Vs + Vist Te” 0, III=(7) 
. ® 
ER 
e(z) = e(z=0) _ { (z' + Ee Se) d2' en Ill=(13) a 
o r 
z err 
8,(z) = —{v, B, — vB) dz', III=(14) 
ki & ® z 
o c 
Vergelijking ke lossen we op voor v $ we vinden dans 
B} 
ls eet 
Vv = ( r áne ) dyste (1) 
We definiëren nu de dimensieloze coördinaat t: t = z/h. Hierin 
is h de schaalhoogte voor de dichtheid in de schijf (zie pag,33 
ze dee = zeer) IIl-(16) 
a GM 
Vergelijking III-(13) wordt: 
t ki B | 
ele) = elt=o) = ( (2tto + prin Lp) ot. (2) 
o 
Voor v, ín vergelijking IIl=(14) kunnen we schrijven (zie pag. 
10 en 31): v (z) 
v_ = & v_(z=0) $ (3) 
1 OE RE v‚(z=0)" 
We substiuëren vergelijking (3) in vergelijking IIl=(1á4) en 
gaan over op de dimensieloze coördinaat t; we vinden zo: DN 
t ART v_(z=0) 
ff hr 
4e) = ( eff ov (er)(til Ara, - B) dt. (4) p 
$ a e° & FE Vg z=0) & z 


Vervolgens definiëren we de volgende constanten: 


1 eg 
e1 = _, C4 = ee, 
RT ATO ee 
v_(z=0) 


hr 
CS = = v‚Cz=0)" C6 = e en C8 = Sr 
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Met deze constanten herschrijven we vergelijkingen (1), (2) en 


(4) en we vinden zoz 


B 2 

v(t) = (C6 + ze) — C8, | (5) 
. 

B‚(t) = (vo Certercsn, _ B_)/C4 dt! en (6) 


ele) =elteo) - ((zergten ver Beer. 
ee 


We kunnen van deze vergelijkingen een eenvoudig blokschema ma= 


kens 


B 
== 2tp - C1 zE Se, 
e(9) 


e= [lers po 


IA 
crjo 





In appendix A ís uiteengezet hoe hiervan een CSMP programma 
gemaakt kan worden, Hoe de diverse constanten berekend moeten 
worden ís uiteengezet op pagina 30 t/m 33, Het uiteindelijk 
resulterende CSMP programma ís op de vogende pagina gegeven, De 
betekenis van de diverse variabelen geven we hieronders 
parameters: M= mj DM = Âj MU = Mons FASTN= Wj Y= r/r, en C 
is gedefiniëerd op pagina 3á, et 
constanten: R/ = ro} RHOP = p(z=0); BZ = B; C3= 1/AN en de 
overige constanten als op de vorige pagina, 
variabelen: TIME = t; VPHI = Vo3 RHO = @; BPHI = B GAMMA = Ys 
LNRHO = Inp; DLNRHO = on en DBPHI = a N 
N.B. om berekeningstechnische redenen wordt.bij de berekening 
van de dichtheid de logarithme gebruikt: op deze manier is ge= 
garandeerd, dat de dichtheid positief blijft, 
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tete ek rk ee 
BEREKESIJG VAS DE “AZINALE GROOTTE VAN HET AZTAIZNAGAT 
HAGNEETVELD IN DE GROBISCHIJF BIJ ES AEETROHENSTER. 


AANMANESS=ALPHIASSCHIJFSODER (SAAKURAASUUYAZY 
„AET be HAMER TS rare HORAE IV DE Meat IE 
EEN FACTOR 0,2 AFGESGAERAS (Gais Ie 


Het AVERGA IGSGEBTES, 
Ton Ysi0h, SIT VELD 
RE DIPIL. CHE FISDES 

EERE 


Ee 
„42 VERDIDERSTELLEN, DAT 
DAT ZICH HITSTREKT VAT 
ZICH GENRAAGT pn Ki 
DAN EEN AOVEAGRINS Va 





DN 
Te a en 
EN be EN 
ei 
De Ì 


TH ELKE RU AOETES DE VOLGENDE PARASETERS HORDEN SEGEVENS 

“HASSA VAN DW HEUTRIHENSTER CM) IN fed ZOISHASAASS 
=AAGSETISCHOIPOJLSOSET CSU) IH 130 GAUS AAK) 

„-FASTJESDPARAHETER CPAST.IJ 

=HASSAOVERDRACHT (22) IN 1517 GRAM/SEC 

„TE ONGERZIEKEN PASITIE CT). TI RO CIS Gha IJ 

=DEGINSAARDE OICHTA0ID CC) Ii DICHT, ALPHASACDEN Se3 
Ii HET V,Be AIERGJDER ZIJN DEZE WAARDE GEGEVEN AUS VOOR 


Kd HH HH He RH HOH OH HM KH HH He aka: 


À Gab IZ) 


Ir 
HERCILES Xe E C&J 
kikkekees tankktkernkt 


KREK KEK KKK 
IGITLAL 
PARAMETER Hale, Hil=3,47, FASTUSI.35, Duster 1250, Csle 
EGTERACTIVE 

Bkk kk Keten kk Eken kde KR KR ek NAR ER KIR Ek kiek eek kk 
* ALS DEZE PARANETERS GEGEVEN ZIJ!, KURKEN SET AEAULD HTERYAU 

& DE CONSTANTEN HORDEN BEREKEND, DIE SE BIJ ON INTEGRATIE 

XK NIODIG HEBBEN, DIT GABENRT AIEROSDERe. 

KRAKEN EERE dr dek ek ee Ad 
Rl TEZKTUKK CAST ed IERCH ATIE TERI STe) 

iOzte JELAKAKRCS O/D EDER A2 IIA CRORIDKECLS SEI EL 

BAse0 2EIIKA/ CROKY)AFF, 

Cist, MIJA 20/1, ZAO) ADAEK(30/ LD, PERIKELEN (IAA) 
C357,25775e2 

CAZ6e SABLEASKAOKNDEN CL /LG ID /DKECET ADRIA) ECS, 116) 


Hit 


% kk KH H& ef 


C3=n3,4846Beh/ AKKI GARNEER CI ATO CI KCROELIARCT ARI Cie SE AGT EELe 5) 


CS53, 2031E26R4/RO/Y 


€ 


Câsi, 78ATELJKSART CARAFE ASTAKY 

Nkk Kekke Sel Bot kk oke ke ee ekke KK OR AK Kk 
_X HJ ALLE PARAMETERS HN CONSTANTE REKESD ZIJ KAS DE id 

* INTEGRATIE HORDEN UITGEVOERD. % 

KRAS ERN KKK KE kk ek le A A kN VARA Ak dk ee ie 

AFHAMIC 


DUIRHOsm 2, KEINE RIZAADITKADRDIJTARHI 
VPIII=ZSARTCCHFCIKAPHETERI ANN) CA 
DIPHISVPHIKCRIAEKTSEBRHTe 2) /04 8 
WI2HOZEITGRL (De, DARO) 

REJZRIOOKEKP CLARE 

APAISINTGRL(Oe, POPHIJ) 

SA AMAZARS(SPII/32) 

TAAMTHAL 

2I4ER DELT20 0991, FISTIN=,, ONTDELSN GOS, DENNAKSD 001 
AZTHON RAS 

FINS V2HI=O, 


PUSH RAI=0, 

NITPJT RA) , BPA , VOET 
Eis KEN 

are 


% HH HH HH od HA HR RHHHR HE HH UH % 
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Appendix Ce. 

In deze appendix geven we het CSMP programma voor de bere= 
kening van de verticale structuur van de accretieschijf in het 
model, waarin het schijfmodel van Ichimaru, 1977, wordt ge 
bruikt, Dit model is beschreven ín de paragrafen IV=1 en 2, De. 
vergelijkingen die we met dit programma willen oplossen, zijn 
dezelfde als in de vorige appendix, De methode van oplossen en 
het programma zijn dan ook geheel analoog. 

Het etoe verschil ín de vergelijkingen ís de radiële druk 
gradiënt in de vergelijking voor v,, Deze vergelijking luidt nu 


Pp 
2 
B } 
cm, Ee ‚B dB 
vs = (5 + me * ame or) SST (1) 


Deze extra term wordt berekend met de methode zoals die op page. 
465 en 46 is uiteengezet, Dit betekent, dat we de berekening 
simultaan op drie verschillende posities le = r= ór, r en 
r. e= rr +óÓr) moeten uitvoeren, waarbij de radiële drukgradiënt 


d 
op deze drie posities gegeven wordt door: 


dB 4B,(r) — 38,(r ) =- B‚(r.) 
B lr) ee, Tr) + + DE, (za). 


d8_ - Br.) -Blr ) 
B Lr) = B rr) + ha a en B en (2b) 


z d 2sr 
dB 3Blr.) + Opr ) - 4B4(r) 
dB z Pe 4 PN b 
Ee CE OS en 


De eerste term aan de rechterkant van deze vergelijkingen is 
onafhankelijk van z en volgt direct uit vergelijking IV=-(19), 
De tweede term aan de rechterkant moet bij elke integratiestap 
opnieuw worden berekend, 

Het CSMP programma, waarmee we de vergelijkingen oplossen, 
bestaat nu uit drie blokken, die elk de berekening voor één van 
de drie posities uitvoeren en vrijwel gelijk zijn aan het 
DYNAMIC deel van het in de vorige appendix besproken CSMP proe= 
gramma, en een vierde deel, dat de relatie tussen deze blokken 
beschrijft en de radiële drukgradiënten op de wijze van verge= 
lijking (2) berekent, 

Ook de gebruikte parameters, constanten en variabelen zijn 
op dezelfde manier gedefiniëerd als in de vorige appendix, Met 
“de achtervoegsels M en P worden steeds de diverse grootheden op 


respectievelijk de posities r en r, aangeduid, 
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De nieuwe constanten en variabelen, die in de vorige ap= 
pendix niet zijn gedefiniëerd en ook niet volgen door de simpe= 
le achtervoegsels M en P,‚ zijns 


parameter: DYREL = 6r/r. 





constanten: OBZDY 


dB_/d(r/r _)s VPHIJ = va(r)3 
en deze beiden met de achtervoegsels M en P, 


variabelen: BOBDY B da/d(r/r) ook met achtervoegsels M en P, 


Andere nieuwe elementen in dít CSMP programma vormen de drie 





statements in elk blok die er voor zorgen, dat als v, in een hal 
positie nul wordt of van teken verandert, Va daar nul blijft, | 
Dit omdat Vp = Q de grens tussen de schijf en de starre magne= 
tosfeer definiëert, zodat voor grotere z we in de magnetosfeer 


zijn, waar v, nul moet zijn, De CSMP statements hiervoor zijns 


& 
(B.V. in de positie r) 


VV=SQRT (C6+(BPHI3#2.+BDBDY )/á ./PI/RHO)-CB 
HH=NOT (VV/VPHIJ) 
VPHI=VVa( 1 „=HH) 


Hierin zijn VV en HH twee variabelen, die we ín deze statements 
slechts als hulp gebruiken om v uit te rekenen. Het NOT staten 
ment is in appendix A gedefiniberd, In de andere posities ge= 
bruiken we dezelfde drie statements, nu met de hulp van VM en 
HM of VP en HP, S 

Op de volgende pagina's is het resulterende CSMP programma 
gegeven, De numerieke waarden van de diverse constanten zijn 
berekend uit de Ichimaru schijfoplossing, dei, vergelijking IVe 


(11), en op de manier zoals op pagina 45 is uiteengezet, 





Br Seker er er EK EE KAR KAK K A9 
% BEREKENING VAi HEI AZINUTRALE MAGUETISCHE VoLu & 


* IN DE GROEISCHIJF BIJ EZ MEUTEONENSTER,. IN DUZE 3Ee * 


* REKENINGEN GEGRUIKEN WE HET HoDub VAL ICHIKARU VOOR « 
* De SCHIJF; IN DAT HODEL IS NANEULIJK OI9K HET EPFFEC= & 
* TIVE GELEIDINGSVERMUGEN BEKEND. & 
* OVERIGE AANGAHESs * 
* „IN DE "Bel." WORDT EVEUALS IN HET LIJDEL VAA G&L + 
* HET VERTICALE VELD EEN FACTIR 0.2 AFSESCHERNC, * 
% * 


* DE TE SPECIFICEREN PARA4ETERS ZIJils HN, D'!, HU, CASTU, & 
* Le DAREL, C 4 
Bits kkk ok ek OE A EE 
INIEIAL 
PARAMETER Hat., Diiet., nU=te., FASIN=0.5, Y=le2e DYRCL=l.E=3, Cale 
INTERACTIVE 
Bkk kkk ek A KA jk 
* BEPALING COUSTANTEN Ii PGSIIIE "Y" * 
Akki keg ek kkk ek bek kkk kk kk 
PI=3.,1416 
RO=1,TEBKHUAR(4/Te D/K 0/7 I/DHEEC2S/ Te) 
RHCO=Le12STELOKNK SSHDUEK A/(ROKY)ARI. SEL 

En ZEZOHKIJ/ROKKZ (119 5/Y*H3e=. 195TEYKRLT3) 
DEZUY=. 2EZOKHU/ROKKZ ACI SBT/YERA, +. 1429/{*4.27) 
C1z4,6OT6E=22/NER I/DijkKe SECROKY) EE, 
C4=3, 3505*5RT(CI) oe | . 
C5=3, 134BE-B/ NERI OSEDMKR O/(RORY) ER I/ (lee YEL. 5) /FAGTE 
C6=3.2031E26*H/RO/YT 
C6=1.78ITELZESORT(Y/RIJKFASTNEY 
VPRIO=SORT(C6+BE*OBEDEY/ 4. /PI/RUII)=C8 
Bkk kk de kk 
X BEPALING COUSTAUTEN IK POSITIE "Zult" p 


eee de eek 








RlOOM=1,125TELOKMKK SSKLMAK A/CROKYECLSDYREL)JKKL 55XC A10 


B4s= 2EL3OHHI/ROFKI, CL. 1I6/CY-YEOYREL) KAJ. 195 TKCY-Y*DYRELD KK. 73) | 
DOZCYEi=, 2E30KHU/ROKK3 ACJ. BET/CY-YKDYREL) KEA. te 1429/(Y=YEDVREL) KK 27) 
CiM=Cik(l.=DYREL)**,9 

C4ii=3, 3565FSARTCCIN)I/ CL. DYREL)KKL.05 

CSHECS/(le=DYRELD EKL ISAAK 5) / Cler CIRCA. -DYREL))KEL,S) 
COM=C6/(1,=*DYREL) 

COH=CEF(L.=DYREL) … | © 
VPHIOCK=SORT(COM+BZNKDBZDYAR(Y=YKDYREL)/ 4e /PI/RHOOA) CBN De 
son kde kk 

« BEPALING CONSTANTEN IN POSITIE "YPLUS" « 

KEKEN AE OE 
RHOOF=1,1257E1OKMER, SSKDMKK, 4/CROKYECI+DYREL) EKL, 55C 

BZPs=, 2E3OFMI/ROKK3, ACL, 196/(Y+YEDYREL) EJee 195 TK(YHYADYREL)K.73) 
DBZDYE=. 2E3OKKU/ROKKIe4C3.BET/CLHYHOYREL) KEG ete 1429/(Y+YADYREL) #27) 
CIP=CIE(1,+DYREL) **,9 

C4AP=3,356S*SIRT(CIPJ/C1.+DYREL) 1.05 

CSP=CS/ (1e +DYREL)AKI35KCH eYRKLL5I/ (le CYR. +DYREL)D EKL, S) 
CEP=CO/(1.+DYREL) | 

CSP=C64(1,+DYREL) 
VPHIOP=SORT(C6P+BZPADBZDYPKCYHY*DYREL)/4. /PI/RHGOP)=CEP 
KEA OE ek 

* DE INTEGRATIE * 


BR EO ok 


DYHEAMIC 

: DE ë 
DLNE Ke 2 TIME/(1,=DYKEL) 2. 1=CAA/4. /PI/RHOMKBPEIKKDEPKIK 
VEESGRT(C6M+(BPHINA#2,+BDBDYH)/4./PI/RKOY) =CBH e 


HEZNOTCVM/VPHIOH) 

VEPHIMSVER(1. =H) 
DEPKIKSVPHINKCHAX(TIHEKCSAKBPHIA=EZM) 
LNEHONZINTSRL(O,., DLNRM) 
RHOKEPHGONKEXPCLNRHOM) 


BEHIMZINISRL(9,. ,DBPHIK) 


arken ode ok A44 
DLNRKU== 2, *TIME=Ci/4, /PI/RHOFEPHIKDDPKI 
VV=SQET(C6+H(BPHIKK2,+5DGDII/4. /PI/RHO)=CB 

HHESOT(VV/VPHIO) 

VPEI=VV&C1 =H) 

DEPHI=SVPHISCAX(TIHEXCS*BPHI=BZ) 

LLRHO=INTSRL(O., OLNRHO) 

REO=RHIOO*EXP(LNEHO) 

BPHI=INTGRL(O., D8PHI) 

GAKHA=BPHI/BZ 
Kekke ek kj ek 
DLIRP== 2, KTIMERCI.+DYREL) A2, 1=CIP/ 4e /PI/RHIPKBEHIPKDLPEIE 
VP=SQRT(CGP+(EPEIP*#2.+EDSDYP)/4./PI/RHCP)= CBP 
HP=ICT(VP/VPHIOP) 

VEHIF=VE*CL =P) 

DEPHIP=VPHIPSCHP*(TINERCSP*GPIIP=BZP) 

LERKOP=INIGRL(O., DLNRP) 

RHOP=KHOOP*EXPCLNRUGP) 

BPHIP=INTGRL(Q., DEPHIP) 

Bkk ek ke eget kk kk fe ek Ek ek kkk jk Kk 
HEM=BZH*DBZDYN 
BDBOYK=CHBY43PHINK(4, XBPHI=3. FBPHIM=EPHIP)/2./DYREL/Y)*(Y=YADYREL) 
BOBCY=(HZ*DEZOY+BPHIK(BPHIP=BPHIMJ/2./DYREL/Y)*Y 
HEP=BZPXDGZOYP 
BUBCYP=(HSP+BPHIP4(3, #BPUIP+BEHIN=4, KBPHIJ/2e/DYREL/Y)*(Y+YEDYREL) 
en 
TERMINAL | 

TINER DELI=C, 001, EINEIK=4., OUIDEL=. 01 

METHOD KKS 

FINISK VPHI=O, | 

OUTPUT BPHI, GAKMA, RIO, VPHI, GPHIK, BPHIP, GZ, ï, Cc 


ELD 





